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This  work-  was  performed  under  US  Army  Natick  Laboratories  Contract 
No.  DAAG17-70-C-0127  during  phe  period  of  1  Apr  70  to  31  Mar  71,  and 
follows  the  work  reported  on  in  report  reference  No.  7*  The  Project  No. 
was  1F162203D195  entitled  "Exploratory  Development  of  Airdrop  Systems", 
ar.d  the  Task  was  No.  13  entitled  "Impact  Phenomena".  Messrs.  Edward  J. 
Giebutowski  and  Marshall  S.  Gustin  of  the  Airdrop  Engineering  Laboratory 
served  as  the  Project  Officers. 

The  effort  is  part  of  a  continui..g  investigation  directed  toward 
obtaining  a  better  understanding  of  the  failure  mechanism  of  energy  dissipater 
materials,  and  the  response  of  airdroppable  supplies  and  equipment  to 
airdrop  imoact  phenomena ;  and  toward  obtaining  improved  ainirop  energy 
dissipater  materials  and  techniques . 

This  report  is  concerned  primarily  with  the  development  of  design 
analyses  of  the  dynamic  response  of  a  complex  structure  to  the  shock  of 
a  vertical  and  planar  airdrop  impact.  Primary  consideration  was  given  to 
the  use  of  the  finite  element  method  of  analysis,  of  a  three  dimensional 
lumped  mass  model  of  the  Army's  M-37  truck. 
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ABSTRACT 


The  principal  purpose  of  this  study  has  been  to 
develop  a  procedure  for  analyzing  the  response  of  complex 
structures  to  impact  and  to  provide  a  computer  code  for 
making  the  necessary  computations.  Attention  has  been 
focused  specifically  on  the  displacements,  velocities  and 
accelerations  produced  at  various  points  in  military  vehicles 
subjected  to  ground  impact  in  airdrop  operations. 

The  vehicle  is  modeled  by  a  lumped  parameter  (spring- 
mass)  system.  Using  the  finite  element  method,  a  set  of 
equations  of  motion  is  formulated  for  this  model.  Then 
these  equations  are  solved  numerically  by  the  Runge-Kutta 
method.  A  model  for  representing  a  specific  vehicle,  namely 
the  M-37  military  truck  cushioned  for  an  airdrop  is  used  to 
illustrate  the  procedure.  Some  of  <the  physical  constants 
for  the  model  are  modified  as  required  to  bring  the  com¬ 
puted  displacements  and  accelerations  'at  various  points  in 
the  model  into  agreement  with  measured  values.,  It  is  found 
that  the  response  of  a  structure  oroperly  cushioned  and 
subjected  to  impact  loading  is  not  sensitive  to.  the  elastic 
properties  of  the  interconnecting  members.  Thus,  the  develop¬ 
ment  of  a  suitable  lumped  parameter  model  of  a  given  vehicle 
is  simplified.  However,  special  attention  should  be  given 
to  the  more  important  components  of  the  vehicle  such  as,  for 
example,  the  engine. 

Experimental  results  show  that  more  than  one  fourth 
of  the  system  energy  is  dissipated  through  the  structural 
damping.  Hence,  damping  must  be  included  in  the  equations 
of  motion. 

Thi  most  essential  factor  affecting  the  dynamic 
response  of  the  system  has  been  found  to  be  the  force  applied 
as  a  result  of  the  impact.  In  the  example  used,  this  force 
is  applied  by  the  cushioning  system. 


INTRODUCTION 

The  finite  element  method  appears  to  be  ideally 
suited  to  the  problem  under  consideration,  namely,  the  com¬ 
putation  of  the  displacements,  velocities  and  accelerations 
at  selected  points  in  a  structure  subjected  to  an  impulsive 
loading.  To  apply  this  method  of  analysis  one  must  first 
prepare  or  select  a  conceptual  representation  of  the  contin¬ 
uous  structure  as  an  assemblage  of  structural  elements  inter¬ 
connected  at  nodal  points.  The  idealized  structure  is  assumed 
to  be  acted  on  by  external  equivalent  forces  and  to  possess 
equivalent  inertia  properties  only  at  the  nodal  points.  Thus 
the  continuous  structure  is  replaced,  for  analytical  purposes, 
by  a  lumped  parameter  system.  The  accuracy  of  the  predicted 
dynamic  response  of  a  structure  will  depend  on  how  well  the 
structure  is  represented  by  Jhe  selected  lumped-parameter 
model.  It  has  been  reported  that  even  for  a  very  simple 
beam  or  uniform  plate  with  boundary  conditions  which  can  be 
exactly  expressed  in  mathematics,  errors  in  the  predicted 
responses  can  easily  be  as  much  as  35  percent.  Although  no 
investigation  of  the  validity  of  lumped  parameter  models  for 
complex  structures,  such  as  vehicles,  has  been  reported  in 
the  literature,  it  may  be  assumed  that  discrepancies  of  even 
more  than  35%  can  be  expected  for  poorly  represented  struc¬ 
tures.  The  general  details  of  th|  finite  element  method  have 
been  discussed  in  the  literature."  However,  the  nature  of 
the  model  which  will  best  represent  a  structure  such  as  a 
vehicle,  with  its  varied  elements,  irregular  geometry  and 
discrete  masses,  is  not  at  all  clear.  The  primary  objective 
in  this  study  will  In  to  select  a  model  and  then  determine  by 
computation  and  by  experiment  how  valid  the  model  is.  It  is 
expected  that  some  suitable  rules  can  be  formulated  regarding 
the  representation  of  vehicles  by  lumped  parameter  models. 

A  model  for  representing  a  specific  vehicle,  namely  the  M-37 
military  truck,  cushioned  for  an  airdrop  is  used  to  illustrate 
the  procedure. 

In  the  procedure  followed  for  this  study,  the  lumped 
parameter  model  for  the  M-37  is  first  developed.  Then  the 
equations  of  motion  of  the  model  are  formulated  following  stan¬ 
dard  finite  element  methods.  These  equations,  with  appropriate 
initial  conditions  are  then  solved  using  the  Runge  -  Kutta 
method.*  Some  important  factors  such  as  structural  elastic  pro¬ 
perties,  damping,  and  impulse  loading  which  affect  the  dynamic 
response  are  investigated.  Finally,  an  experimental  prograr 
of  actual  truck  drop  tests  is  carried  out,  and  the  results  are 
compared  with  computed  results.  In  the  analysis  of  the  i athe- 
matical  model,  the  concept  of  linear  transformation  is  exten¬ 
sively  used.  Linear  transformation  techniques  streamline  and 
simplify  considerably  the  procedures  involved  in  the  analysis. 

It  should  be  mentioned  here  that  the  computer  program  devel¬ 
oped  for  the  vehicle  is  also  applicable  to  other  complex  struc¬ 
tures  . 


*  Algorithms  for  this  method  are  available  in  most  computa¬ 
tional  facilities. 
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MATHEMATICAL  MODEL 

A  vehicle  such  as  the  to-37  truck  may  be  represented 
by  the  model  shown  in  Fig.  2.1.  In  this  model  the  engine, 
transmission,  transfer  case,  differentials  and  wheels  ars 
treated  as  discrete  masses.  The  mass  of  the  winch  is  assumed 
to  be  distributed  uniformly  along  the  two  main  longitudinal 
and  the  remaining  transverse  members  of  the  truck  frame.  The 
adoption  of  this  model  is,  however,  quite  arbitrary.  Many 
other  arrangements  of  masses  wculd,  no  doubt,  be  equally  accept¬ 
able  . 

When  a  structure  such  as  this  vehicle  is  to  be  intention¬ 
ally  subjected  to  an  impact,  as  in  an  airdrop,  cushioning  is  pro¬ 
vided  to  reduce  the  severity  o.‘  the  shock  produced  by  the  impact. 
Usually  all  of  the  discrete  masses  shown  in  Fig.  2.1  would  be 
cushioned  independently,  if  possible.  For  the  M-?7  truck  the 
engine  is  not  cushioned  independently,  partly  because  i;  is  shock 
mounted  on  the  frame,  and  partly  because  of  geometrical  and  struc 
tural  problems.  The  engine  i a  supported  on  rubber  cushions,  or 
shock  mounts  at  three  points,  one  in  front  and  two  in  the  rear. 
The  action  of  these  mounts  can  be  represented  by  the  spring- 
damper  system  shown  in  Fig.  2.2.  The  stiffness  and  the  damping 
capacity  of  a  mount  depends  on  many  factors  such  as  the  hardness 
of  the  rubber,  the  shape,  and  the  age.  Rather  than  try  to  deter 
mine  a  precise  set  of  values  for  the  M-37  mounts,  values  of 
20,000  lb/in.  for  k  and  60%  of  critical  for  the  damping  were  arbi 
trarily  assumed.  Later  these  values  were  varied  to  improve  the 
"fit"  between  experimental  observations  and  computed  results. 

The  transmission  is  actually  attached  to  the  engine 
but  since  it  can  be  cushione 1  independently,  it  has  been 
assumed  to  be  a  mass  which  is  attached  to  the  engine  by  a 
very  stiff  element. 

The  transfer  case  is  supported  at  four  points  on  two 
central  cross  frame  members  which  act  as  spring  supports  for 
it. 

Wheels  and  differentials  are  connected  by  the  axles 
which  in  this  analysis  are  assumed  to  be  rigid,  massless  rods 
attached  to  leaf  springs.  Tires  absorb  considerable  e  lergy 
and  this  energy  is  given  back  in  rebound.  Also  measured 
relative  displacements  between  the  rear  axle  and  truck  frame 
indicate  that  little  energy  is  dissipated  through  the  5our 
shock  absorbers  associated  with  the  wheels  during  the  Impact. 
Consequently  the  shock  absorbers  are  neglected  and  the  vehi¬ 
cle  is  supported  by  the  leaf  springs  in  Fig.  2.3. 

The  spring  constants  K  and  damping  factor  C  for  tne 
tire  in  Fig.  2.3  must  be  determined  experimentally,  or  be 
estimated  using  whatever  guidance  is  available.  These  quan¬ 
tities  are  initially  assumed  to  be  7,000  lb/in.  and  2C%  res¬ 
pectively  . 

The  winch  and  the  other  distributed  masses  of  the 
truck  are  cushioned  with  the  two  main  longitudinal  frame 
members.  in  the  model  shown  in  Fig.  2.1,  all  individual 
components  are  interconnected  at  the  nodes  numbered  from 
1  to  39.  In  the  present  study  only  a  vertical,  planar  im- 


2 


Fig.  2.1  Lumped-Parameter  Model  -  M-37  Truck. 
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Fig.  2.3  Suspension  System  of  Vehicle 


Fig.  2.4  Element  Jisolacem^nts 


pact  is  considered.  Consequently  each  nodal  point  in  the 
model  is  to  have  onlv  3  degrees  of  freedom,  consisting  of 
one  vertical  translation  and  two  nlanar  rotations.  Three 
degrees  of  freedom  are  associated  with  each  end  of  all  members 
as  shown  in  Fig.  2J*. 


ANALYSIS  OF  THE  MODEL 

For  the  model  shown  in  ?-'ig.  2.1  with  assumed  damning 
subject  to  imoact  loading,  the  eouation  of  motion  can  be 
written  as: 


1  +  £  i  +  %  2.  "  £(t) 


(3.1) 


where  M,  0  and  £  are  square  matrices  of  inertia,  damning 
and  stiffness  respectively  and  F ( t )  is  the  column  matrix 
of  cushion  forces.  The  generalized  displacement  matrix  q_  is 
numbered  in  the  senuence  according  to  the  joint  numbers.  For 
example,  coordinates  .  _  1  and  q,.  are  associated  with 

the  translator  and  tworotational  motions  about  common  datum 
axes,  at  the  i^h  nodal  joint  of  the  model.  For  example,  at 
Joint  5,  the  coordinates  are  ^^-2  I  q ^ ^  q  or  (]n  o,  ^ 

and  q  ^  . 

The  stiffness  matrix  K  of  the  comnlete^element  assembly 
Is  obtained  by  the  direct  stiffness  method.  This  method 
consists  of  first  deriving  the  individual  element  stiffness  as 
in  element  coordinates,  followed  by  a  coordinate  transforma¬ 
tion  and  the  subsequent  superposition  of  each  element  stiff¬ 
ness  so  that  the  translational  and  rotational  degrees  of  free¬ 
dom  of  all  elements  which  share  a  common  nodal  point  are  ex¬ 
pressed  in  the  same  coordinates.  The  superposition  of  each 
transformed  stiffness  is  accomplished  by  adding  its  individual 
terms  into  the  complete  stiffness  matrix  according  to  the 
element  nodal  point  numbers.  The  same  method  can  be  emploved 
to  obtain  the  mass  matrix  H  and  the  force  vector  F.  The 
derivations  of  mass  matrix"M,  stiffness  matrix  K,  force  vector 
F  and  damping  matrix  g  are  discussed  in  detail  as  follows: 

1.  Element  Displacement  Functions 

To  determine  the  elastic  and  inertia  properties  of  a 
structural  element  the  strain  and  kinetic  energies  of  the 
actual  continuous  element  are  equated  to  the  corresponding 
quantities  for  the  equivalent  discrete  model.  Consider  now 
a  uniform  structural  element  in  the  horizontal  plane  as  shown 
in  Fig.  3.1-  The  common  datum  (X,Y,Z)  is  established  for  all 
structural  elements  so  that  all  displacements  and  corresponding 
forces  will  be  referred  to  this  common  coordinate  system.  The 
origin  of  element  coordinates  (x,;/,z)  is  located  at  node  1 
with  the  ox  axis  taken  along  the  length  of  element  and  with  the 
oy  and  oz  axes  as  the  orincinal  axes  of  the  element  cross 
section. 

The  column  matrix  U  for  this  element,  as  mentioned  before, 
consists  of  six  displacements,  two  vertical  (Z  direction)  deflec¬ 
tions  U1  and  Ujj  and  four  rotations,  >  ^5  j  an^  dg- 
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The  displacement  furction  w^  may  be  assumed  as 

wz  Cl  +  C2X  +  C3X  +  C‘ix^  +  C5y  +  C6xy 

where  C,  to  Cg  are  constants  to  be  determined  by  boundary 
conditions  at  Doth  ends  of  the  element. 

At  node  1:  (x  =  0,  y  =  0) 

w  =  U  =  Cn 
z  1  1 

3w 

I3T  "  "i;3  =  C2 


_ ~  =  +  it  =  r 

9v  U2 

At  node  2:  (x  =  L,  y  =  0) 


w  =  U,  =  C.  +  C  „L  +  C-l/  +  C  i, 
z  4  1  2  3  4 

9wz  2 

I3T  "  -  U6  "  C2  +  2CJL  +  3C4L 


U5  =  C5  +  C6L 


Solving  for  the  C's, 

C1  =  U1 
C2  "  -  U3 


c  =  — - 

3  L2 


[  -  3’^  +  2LU3  +  3U4  +  lu6  ] 


Ch  =  ^  [  2UX  -  LU3  -  2U4  -  LUg  ] 

C5  "  T;2 


n 


[  -  3-  Uc  ] 


Define  the  nondimensional  parameters 


x 

r  =  r 


■■ l 


Thus  the  displacement  function  w  becomes 

w  =  (1  -  3r2  +  2r3  )U,  +  Ls  d  -  r)U0 
z  1  z 

+  L(-  r  +  2r2  -  r3)U3  +  (3r?  -  2r3)U4 

+  LrsUc  +  L(r2  -  r3)U, 

5  6 

The  displacement  function  w^  caused  by  twisting  is 


(3.7) 


(3.8) 


wr  =  -  Ltd  -  r)_U2  -  LrtU5  (3.9) 

This  is  a  direct  geometric  relationship,  .it  can  be  obtained 

by  assuming  =  c  .  c  x  +  r  x2  +  C,x3  +  CRz  +  Crxz. 

y  1  2  3  4  5  6 

The  displacement  function  w  is 
9  w  9  w  X 

w  =  -  z-J.  -  y-4  (3.10) 

X  9x  1  9x 

Then  by  combining  3.3,  3.9,  and  3.10  obtain 
w  =  6t(r  -r2)U  +  Lt(l  -  4r  +  3r2)U 

A  J.  J 

+  6t(-  r  +  r2)U4  5-  Lt  (-  2r  +  3r2)L’6  (3.11) 

Thus  the  displacement  functions  w  can  be  expressed  in  terms 
of  the  discrete  displacements  U  as 

w  =  a  U  (3.12) 

where  a  is  the  3x6  matrix  whose  transpose  is  given  by: 


T 

a  = 


6t(r  -  r  ) 

0 

Lt(l  -  4r  +  3r2) 
6t(-  r  +  r2) 

0 

Lt ( -  2r  +  3r2 ) 


0 

-  Lt(l  -r) 
0 
0 

-  Lrt 
0 


2  3 

1  -  3r  +  2r 


-  Ls  (1  -  r) 

L (-  r  +  2r2  -  r3) 

2  3 

3r  -  2r 


-  Lrs 

T  ,  2  3 . 

L  ( r  -  r  ) 


(3 


2.  Element  Mass  Matrix 

To  determine  the  element  mass  matrix  begin  with  the 
element  displacement  functions  w  in  terms  of  discrete  dis¬ 
placement  U  as  indicated  by  E  q.  3.12. 

The  transpose  of  w  is: 


T  rlT  T 
w  =  U  a 


(3.14) 


13) 
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T  T 

=  6U  v  p  a  a  dV) U 
>  V - 

where  the  integration  is  performed  over  the  whole  volume  of 
the  element. 

Hence  the  element  mass  matrix  is  defined  as: 

T 

m  =  p  a  a  dV  (3.18) 

-  J  v - 

Substituting  the  matrix  a  [Eq.  (3.13))  into  Eq.  (3.13)  and 
performing  the  integrations ,  the  resulting  6x6  mass 
matri.v ,  Eq.  (3.19)  is  obtained.  In  the  mass  matrix,  terms 
with  the  moment  of  inertia  I  represent  rotatory  inertia  and 
terms  with  the  polar  moment  Xf  inertia  I  represent  the  tor¬ 
sional  inertia  of  the  element.  From  the  results  of  numerical 
computation,  the  effects  of  rotatory  ana  torsional  inertia  of 
the  element  during  impact  are  considered  negligible. 

3.  Element  Stiffness  Matrix 

The  bending  strain  in  the  xz  plane  of  the  element 
shown  in  Fig.  (3.1)  is: 

2 

3  w  ,  1 

e  =  - 1  =  -0(-  6  +  12r) U,  +  =-(4  “  6r)U, 

D  2  _  z  1  L  j 

3X  L 

+  — „ ( 6  -  12r)U.  +  ~(2  -  6r ) Ur  (3.20) 

2  4  L  b 

Li 

The  twist  strain  is: 

et  =  E<-  U2  T  U5>  (3'21) 

The  total  strain  e  in  matrix  form  is: 

e  =  b  U  (3.22) 


where  e 


8 


symmetric 


156  +  504— _ 

AL2 

I 

0  140-i 

A 

1  ,  1 
-22L  -  42r-^  0  4L  +  56-4 

AL  A 


m  - 


pAL 

T20 


I 

54  -  504-^ 

AL2 


-13L 


+  4  2^ 

AL 


I 

156  +  504-X- 
AL2 


0 


I 

140- 


-13L 


42Sl 


-3L' 


-  14-4 

A 


22L  +  42^ 


4L* 


(3.19) 


+  56 
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>1* 


1  +  2r)  0  £(4  -St)  — ~  ( 6  -  12r )  0  i(2  -  6r) 

l2  L  ir  L 


(3.23) 

The  bending  moment  ip: 

S.  =  El  e, 
b  y  b 

where  E  is  elasticity  modulus.  The  twisting  moment  is: 

St  =  GJet 

where  G  is  the  shear  modulus  and  J  is  the  torsional  constant 
for  the  element. 

Thus  the  stress-strain  relationship  in  matrix  form  is: 

S  =  A  e  (3.24a) 


Where : 


(3.24b) 


(3.24c) 


From  Eq.  (3.22),  we  have  the  virtual  strain 
Se  =  b<5U , 

and  the  transpose  of  virtual  strain 

T  T  T 

6e  =  6U  b  . 


Thus  the  virtual  strain  energy  expression  is: 

L  T 

<5V  ,  .  =  Se  s  dx 


strain 


(3.25) 


(3.26) 


6UTbTA  b  (J  dx 


T  ^  T 

=  6U  (  dA  b  dx)U 


(3.27) 
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The  element  stiffness  Is  then  defined  as: 


k 


yTA  t  dx 


(3.28) 


Substituting  the  matrix  [Eq.  (3.23)]  and  4  [Eq.  (3.24c)]  into 
Eq.'  (3-28)  and  performing  the  integrations,  the  6x6  stiffness 
matrix  k  is  obtained  as  follows: 


2 

GJL 

El 

y 


o  6l 


symmetric 


12 


(3.29) 


GJL2 

El 

y 

0 


0 


0 

6l 


4L 


2 


i-  J 

4.  Element  End  Forces 

In  general,  the  virtual  work  done  by  the  external  forces 
n(x,t)  acting  along  the  element  at  any  instant  in  time  through 
virtual  displacements  6w  is: 


^external 


jf 

|  6ir 

Js 


T 

6w-£  ds 


T 

|  £  ds 


6UT  t  |T£  ds 

*  C 


Thus  the  eouivalent  concentrated  forces  f  at  the  element  ends 
due  to  a  distributed  loading  £  is: 


(3.30) 
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However,  the  forces  f  in  Eq.  (3.30)  are  the  fixed  end  forces 
(with  the  signs  reversed)  of  the  element  corresponding  to 
the  distributed  forces  £  along  the  element.  Two  cases  of 
cushioning  forces  applied  on  the  element  as  shown  in  Fig.  3.2 
are  considered  in  the  analysis. 


The  fixed  end  forces  for  case  1  are: 


where : 


f  = 


i  PLX  +  Ml 


M- 


L 

0 

-M, 


PLj  +  M2 


-  M, 


M,  = 


M-  = 


PL1L2 


PL1L2 


L 

0 

M. 


(3.31) 


For  case  2, 


f  = 


where : 


[?<L 


L]_  +  L2)  (L  -  Lx  -  L2)  +  M3  -  M4]/L 
0 


[^(L  +  L1  -  L2)  (L  -  L1  -  L0)  +  M4  -  M3]/L 

0 


+M, 


M-.  = 


— ,  [  (L  -  L1  )  3  (L  +  3L,  ) 
12L2  1 


(3.32; 


L2  (4L  -  3L2) ] 


M 


— ,  [ (L  -  L-) 3  (L  +  3L0 ) 


L?  (4L 


5.  Linear  Transformation 

If  the  set  of  coordinates  q  with  n  degrees  of  freedom 
in  the  equation  of  motion  Eq.  (3.1)  is  a  linear  combination  of 
a  different  set  of  coordinates  u  with  m  degrees  of  freedom, 

2  =  B  u  (3.33) 

then  the  mass,  damping,  stiffness  and  force  matrices  in  u 
coordinates  can  be  calculated  by  using  the  concepts  of  energy 
and  virtual  work.  The  elements  of  the  n  x  m  matrix  B  [Eq.  (3.33)] 
are  constants. 

In  the  q  coordinate  system,  the  kinetic  and  potential  energy 
expressions  have  the  matrix  forms: 

T  =  j  4TM  2.  (3.34) 

v  =  \  aT|  2  (3.35) 


where  M  and  K  are  mass  and  stiffness  matrices  in  q  coordinates 
re  specTTi  ve  lyT 
From  Eq.  (3.33)  , 

£  =  B  u  (3.36) 


end  the  transposed  matrices 

T  T_T 
2.  ~  2  ® 


.T  . TnT 

2  =  2  B 


(3.37) 

(3.38) 


Introducing  the  linear  transformation  Eqs.  (3.33),  (3.36),  (3.37), 
and  (3.38)  into  the  kinetic  and  potential  energy  expression  Eqs. 
(3.34)  and  (3.35)  results  in: 


T  = 


1  aTm  • 

2  4  M  2 

1  .  T 

U  B  M  B  U 

1  -T  • 

2  -  =  - 


(3.39) 


V 


1 

2 

1 

2 

1 

2 


T 

2  |  2 

T  T 

u  B  K  B  u 

t" 

u  k  u 


where.  T 

m  =  B  M  B 

T 

k  =  B  K  B 


(3.40) 


(3.41) 

(3.42) 
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are  the  corresponding  mass  and  stiffness  matrices  in  the  u 
coordinate  system.  Since  M  and  K  are  symmetric,  it  follows 
that  m  and  k  are  symmetricT  The-"virtual  work  done  by  the 
external  forces  F  in  q  coordinates  acting  through  the  virtual 
displacement  6c[  Ts: 

6W  =  6£TF  (3.43) 


Since  the  virtual  displacements  6a  are  related  to  the  virtual 
displacements  6u  by: 


6q  -  B5u 

(3.44) 

and : 

_4.  —  — 

T  T  T 

62  =  «u  B 

(3.45) 

thus 

T  T  T  T 

5W  =  62  F  =  5u  B  F  =  6u  f 

(3.46) 

where 

T 

f  =  B  F 

(3.47) 

are  the  external  applied  forces  corresponding  to  the  u 
coordinates . 

Consider  now  the  viscous  damping  force  D.  These  internal  forces 
may  be  expressed  as: 

D  =  -  C  2.  (3.48) 


where  C  is  the  damping  matrix  in  q  coordinates.  The  virtual 
work  done  by  the  damping  forces  D  is: 

5W  =  02Tn>  =  -62TC  2  =  -SuT3TC  B  u  =  -6uTc  u  (3.49) 


v*  i. in—  j.  ^  rp 

c  =  B  C  B  (3.50) 

is  the  damping  matrix  in  u  coordinates . 

In  summary,  the  mass  matrix  M,  stiffness  matrix  K, 
damping  matrix  C  and  force  vector  F  Tn  the  q  coordinate  'sys¬ 
tem  under  the  lTnear  transformation  2  =  B  u  are  transformed 
into  the  element  mass  matrix  m,  the  stiffness  matrix  k,  the 
damping  matrix  c,  and  the  force  vector  f_  in  the  new  u— coor¬ 
dinate  system  according  to  the  following  transformations: 


T 

B  M 

B 

(3.51a) 

T 

B  K 

B 

(3.51b) 

T 

B  C 

B 

(3.51c) 

T 

B  F 

(3 . 51d) 

6.  Transformation  of  Element  Coordinates  to  Datum  Coordinates 
Since  the  element  mass  matrix  m,  stiffness  matrix  k, 
and  end  force  vector  f_  are  initially  calculated  in  local  ele¬ 
ment  coordinates,  suitably  oriented  to  minimize  the  computing 
effort,  it  is  necessary  to  introduce  transformation  matrices 
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u 


changing  the  frame  of  reference  from  a  local  to  a  datum 
coordinate  system.  Consider  again  the  element  shown  in 
Fig.  3.1.  to  Ug  are  displacements  in  the  direction  of 

local  coordinates  x,y,z  and  to  Qg  are  displacements  in 
the  directions  of  datum  coordinates  XYZ.  0  is  the  angle 
between  ':he  r  and  X  axes.  The  displacements  U  can  be  related 
to  displacements  Q  as: 

U1  =  Q1 

U2  =  Q 2COS0  +  QgSi113 
=  -Q2sin0  +  Q^cosO 

U4  =  Q4 

Uc  =  Qccos0  +  Q,sin0 
b  5  6 


U,  =  -Qcsin0  +  Q,cos0 
b  d  b 


or  in  matrix  form 

U  =  B  Q 

(3.52) 

where 

1 

0 

r\ 

u 

0 

cos 3 

sin© 

B  = 

0 

-sin0 

COS0 

1 

0 

0 

0 

0 

CO  90 

sin0 

-sirD 

COS0 

(3.53) 

The  element  mass  matrix  m* ,  stiffness  matrix  k  * ,  and  the  force 
vector  f_*  in  q-  displacentents  can  be  obtained~t>y  using  the 
transformations  of  Eqs.  (3.51a),  (3.51b),  and  (3.51d)  as: 

T 

m*  =  B  m  B 

k*  =  BTk  B_  (3.54) 

T 

f*  =  B  f 

7.  System  of  Assembled  Structure 

Since  all  element  stiffness  matrices,  mass  matrices 
and  force  vectors  are  now  referred  to  the  common  datum,  the 
stiffness  matrix,  mass  matrix  and  force  vector  of  the  complete 
element  assemblage  can  be  obtained  by  the  direct  stiffness 
method  as  mentioned  previously.  The  concentrated  masses  at 
nodal  joints  such  as  transmission,  differentials  and  wheels 
are  simply  added  to  the  corresponding  diagonal  terms  in  the 
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assembled  mass  matrix.  For  example,  the  mass  of  the  wheel  at 
node  no.  24,*  which  is  associated  with  ■‘■he  coordinate  q,  can 
be  directly  added  to  the  diagonal  element  m.^Q ,  ^  of  the  n.u.st 
matrix  of  the  system.  The  concentrated  cushion  for ces  applied 
,  nodal  points  are  also  simply  added  to  the  corresponding 
terms  in  the  system  force  vector.  For  example,  the  cushion 
force  applied  at  the  transmission  can  be  added  to  the  element 
f g 2  of  the  system  force  vector  F.  "owever,  the  discrete  masses 
such,  as  the  engine  and  transfer  case  which  connect  several 
nodal  points  can  not  be  superimposed  on  the  system  directly. 

A  special  study  of  these  masses  is  necessary. 


8.  Engine  and  Transfer  Case 

Since  the  engine,  transfer  case  and  similar  masses 
are  all  supported  at  several  nodal  points,  special  considera¬ 
tion  of  these  parts  is  required.  Consider  now  a  discrete  mass 
supported  by  n-springs  with  coordinates  q^....q  .  Since  three 
points  define  a  plane,  q^,  and  q^,  as  generaSized  coordinates, 
sufficiently  define  the  motions  of  mass.  Let  u-^,  u„,  and  u~  be 
the  rotations  about  the  principal  axes  of  mass  and  the  vertical 
translation  respectively.  The  relationships  between  coordi¬ 
nates  q's  and  u's  can  be  formed  in  the  following  way: 

Fig.  3.3  shows  the  system  under  consideration: 


al'  a2  =  T^e  uni*-  vectors  along  two  prin¬ 
cipal  directions  of  the  mass. 

** 

l,  g,  k  =  unit  vectors  along  datum  coor¬ 
dinates  X,  Y,  Z. 

Xi'  Yi'  ^i  =  coordinates  of  supports. 

cT  ,  c^,  =  position,  vectors  of  points  where 

the  springs  attach  to  the  mass,  re¬ 
ferred  to  the  center  of  mass. 

e^  =  vectors  from  support  point  q-^  to  q^ . 

For  small  rotations, 

q^  =  u^  (a-,  x  cL  •  it)  +  02(^2  x  ^  •  k)  +  ug 

i  =  1,2,3 

which  in  matrix  form  is 


*1 

x  ^  it  a2  x  5^  •  it  1 

\1 

! 

^2 

= 

al  x  a2  ’  ^  ^2  x  ^2  *  it  1 

U2 

.  q3  . 

al  x  ^3  ’  ^  a 2  x  ^ 3  *  ^  1 

.  U3  . 

or  2.  =  AU 


(3.55) 


*  The  q  coordinates  at  node  are  q-jQi  translation  in  the  Z 
direction,  q.^,  rotation  about  the  x-axis  and  q.^,  rotation 
about  the  y-axis. 
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Thus 


where 


u  =  A  ^C[  =  B 
B  =  A_1 


(3.56a) 


(3.56b) 


If  m-^,  it^  and  m,  are  the  mass  moments  of  inertia  and  the  mass, 
then  the  matrix  of  the  discrete  masses  m  the  u  coordinate  sys¬ 
tem  can  be  represented  as: 


According  to  the  linear  transformation  (3.56a),  the 
mass  matrix  corresponding  to  the  q  coordinate  system  is: 


=  B  m  B 


If  there  are  cushioning  forces  f1 ,  f„ 
in  u^,  u^,  and  u^  directions,  then  th 


and  directions  are: 


(3.57) 

and  f^  directly  applied 
e  forces  in  the  q-,  ,  q~. 


f*  =  Bxt 


(3.58) 


Now  the  mass  matrix  m*  can  be  added  to  the  mass  matrix  of 
the  complete  element- assemblage  and  the  forces  f*  can  be 
added  to  the  system  forces.  It  must  be  noted  that  before 
taking  this  step,  the  constraint  coordinates  q^  through  qn 
should  be  eliminated  from  the  system.  To  do  this,  first 
find  a  transformation  matrix  B  such  that: 


(3.591 


and  then  by  linear  transformation  get  the  modified  mass, 
stiffness  and  force  matrix  of  the  whole  system. 

The  transformation  matrix  B  in  Eq.  (3.59)  is  formu¬ 
lated  by  finding  all  the  expression's  cf  q^  through  qn  in 
terms  of  q^,  < and  q^.  Consider  now  the  displacement  q^ 
as  indicated  in  Fig.  3.3.  The  vectors  from  support  1  to 
2,  3  and  i  can  be  expressed  as: 
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L  *g 


e2  =  (x2  -  xx)  1  +  (y  -  y1)j  +  (cr?  -  q.^  £ 
g 3  =  (x  -  x)l  +  (y3  -  Y1)j  +  (q3  ~  q-^K 


(3.60) 


(3.61a) 


(3.61b) 


x2  -  xi  y2  -  yl 

xt  -  x1  yi  -  y1  (3.61c) 

e3  =  x2  '  xi  y 2  “  yl 

*3  -  x1  y3  -  y1  ( 3 . 61d) 

Thus  all  displacements  through  q  can  be  expressed 
in  terms  of  q^,  q~t  and  q3  simply  by  replacing  the  index  i  in 
the  determinants  and  E2  by  the  numbers  4  to  n. 

9.  Static  Condensation 

If  an  element  is  not  rigidly  connected  to  another 
element,  for  example,  a  hinged  connection,  the  element  stiff¬ 
ness  matrix,  mass  matrix  and  fixed  end  force  vector  must  be 
modified.  If  some  coordinates  included  in  the  static  analysis 
are  excluded  from  the  dynamic  analysis,  or  some  coordinates 
with  zero  or  very  small  mass  must  be  removed  from  the  equa¬ 
tions  of  motion  to  avoid  unreasonable  numerical  results  (infi¬ 
nite  values  in  computed  accelerations) ,  then  the  system  mass 
matrix,  stiffness  matrix  and  system  force  vector  must  all  be 
modified.  All  the  modifications  can  be  achieved  by  a  static 
condensation  technique.  By  using  this  cechnique  a  trans¬ 

formation  matrix  similar  to  the  matrix  B  in  Eq.  (3.33)  can 
be  formulated.  Then  the  modified  matrices  can  be  written  in 
the  format  of  Eq.  (3.51). 

The  first  step  is  to  partition  the  mass  matrix  M, 
stiffness  matrix  K,  force  matrix  F  and  displacement  matrix 
U  into : 
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1 


(3.62a) 


M 

=22 


=22 


(3.62b) 


(3.62c) 


(3 . 6  2d ) 


The  column  matrix  0^  refers  to  all  the  displacements 
we  wish  to  retain,  while  U2  aenotes  all  the  remaining  displace¬ 
ments  which  will  not  be  employed  in  formulating  the  new  equi¬ 
valent  matrices.  The  displacements  U2  may  be  determined  from 
the  static  equilibrium  equation  F  =  K^U  by  assuming  that  the 
forces  F2  corresponding  to  the  displacements  U_2  are  axl  equal 
to  ?,ero.  Hence: 


£2  =  -S22  £21  -1 


thus 


-522  £21 


U]_  =  U]_ 


(3.63) 


(3.64! 


where  I  is  a  unit  matrix  and 


3*  =1  T,-l 


=22  =21 


(3.64a) 


Thus  the  modified  mass  matrix  M* ,  stiffness  matrix  K* ,  and 
force  matrix  F*  for  displacements  U  are  obtained  from  Eq. 
(3.51)  as:  1  T 

M*  =  B*  M  B*  ,  K*  =  B*'K  B*  ,  F*  =  B*  F  (3.65] 


Eq.  (3.65)  may  be  obtained  directly  from  the  equations 
of  motion  as  follows: 


=11-1  +  =12-2  +  =11-1  +  =12—2  -1 


(3.66a) 


=21-1  +  =22-2  +  -21-1  +  -22-2  --2 


(3.66b) 
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Eq.  (3.66b)  can  be  rewriv'j:;:.  a-: 

"2  =  -|i2l2l“l  +  l22(  £2  -  ^2liii  -  522“2  >  <3'66=> 

If  we  assume  F„  and  U  to  be  linear  in  a  small  time  interv.  1, 
then  the  second  time  derivative  of  Eq.  (3.66c)  is 


-2  =  22— ?1— 1  ( 3 . 66d) 

Substituting  Eqs.  (3.66c)  and  (3,66d)  into  Eq.  (3.66a) 
the  mass  matrix  M* ,  stiffness  matrix  K*  and  force  matrix  F*  are 
obtained  as  shown  in  Eq.  (3.65).  These  are  the  matrices  that 
will  now  be  used  in  Eq.  3.1  for  actual  computations. 

10.  Structural  Damping 

Until  now  the  discussion  of  the  damping  matrix  in 
the  equations  of  motion  Eq.  (3.1)  has  been  intentionally 
avoided.  However,  earlier  experimental  studies  of  the  M-37 
truck  indicate  that  an  appreciable  amount  of  energy  may  be 
dissipated  through  internal  friction  within  the  truck  body 
or  at  joints  between  frame  elements.  Structural  damping 
forces  D  are  proportional  in  magnitude  to  the  internal  elas¬ 
tic  forces  and  opposite  in  direction  to  the  velocity  of  the 
system.  The  force  may  be  represented  as: 


D  =  i  bq  K  2  (3.67) 

where  b-q  is  a  constant  and  i  is  the  imaginary  number. 

This  expression  for  the  damping  force  is  not  readily 
amenable  to  structural  analysis.  Thus  the  concept  of  viscous 
damping  is  used  in  the  analysis  as  indicated  by  the  damping 
matrix  C  in  Eq.  (3.1).  Furthermore,  this  damping  matrix  may 
be  assumed  to  be  proportional  to  the  mass  and  stiffness  matrices 
as : 

C  =  a  M  +  b  K  (3.68) 

where  a  and  b  (b  is  not  the  same  as  b-q)  are  constants  and  can 
be  determined  so  as  to  give  reasonable  damping  in  the  system. 

If  Eq.  (3.68)  is  substituted  into  Eq.  (3.1),  the  modified 
differential  equation 

M  2.  +  (a  M  +  b  K)  2  +  K  2  =  F(t)  (3.69) 

is  obtained.  — 

Since  both  M  and  K  have  non  zero  off-diagonal  coefficients, 
these  differential"  equations  are  coupled.  By  using  the  modal 
matrix,  obtained  from  the  arrangement  of  normal  modes,  as  a 
transformation  matrix,  the  equations  of  motion  can  be  reduced 
to  a  set  of  uncoupled  equations  of  the  type: 


m  q  + 
r  r 


^am  + 
r 


bk 


q  + 


k  q 
r^r 


=  fr(t) 


(3.79) 


where  q  is  the  normal  coordinate  in  rth  mode  of  the  system 
and  m  and  k  are  the  corresponding  mass  and  stiffness  in 
that  mode  respectively. 

Since  this  equation  is  the  same  form  as  the  single 
degree  of  freedom  equation,  the  critical  damping  in  the  rth 
mode  is: 


(am  +  bk  ) 
r  r  cr 


=  2m  u 

r  r 


(3.80) 


where  to' 


=  kr/mr  is  the  natural  frequency  of  rth  mode. 
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(  Pi  N 
\  J  •  >  J  ) 


Dividing 


bo  oh 
(a 


sides 
+  bu>r2 


of  Eq. 

), 


cr 


(3.80)  bv  Mr, 

2ii) 

r 


or 


air 

2 


where  S  is  the  ratio  of  actual  to  critical  damning. 

Fo?  given  values  of  a  and  b  the  frenuencv  s  which 
yields  a  minimum  value  for  the  damning  ratio  S  can  be  found 
by  differentiating  En .  ( 3 -  8 1 )  with  resnect  to  w  and  then 
setting  the  derivative  equal  to  zero.  Thus 


and 


(3.82a) 

(3.32b) 


If  S  and  a)  are  given,  the  damning  coefficients  a  and  b 
are  calculated  from  the  following  eouations: 


and 


a  =  S  a) 
b  =  S/u) 


(3.83a) 


(3.83b) 


Therefore,  if  the  significant  frequency  range  is 
for  a  structure  and  the  equivalent  modal  damning 
the  constants  a  and  b  can  be  calculated  from  Eos. 
However,  in  the  numerical  analysis  of  this  study, 
of  the  structure  is  not  calculated.  The  value  of 
mated  as: 

1 


i 


established 
is  selected, 
(3.83)  . 

the  frequency 
id  is  esti- 


(3.84) 


where  m.  and  k.  are  the  diagonal  elements  in  the  n  x  n  mass 
matrix  and  stirfness  matrix. 

It  appears  that  for  the  1-37  truck  model  a  reasonable 
value  of  damping  ratio  S  corresponding  to  the  frequency  u  in 
Eci .  (3.84)  is  around  0.01.  This  number  was  estimated  after  a 
comparison  between  some  computed  and  some  experimental  results 
was  made. 


NUMERICAL  SOLUTION 

1.  Runge-Kutta  Method 

The  equations  of  motion  Ea .  (3.69)  can  be  solved 
numerically  by  using  the  Runge-Kutta  method.  However,  Ea . 
(3.69)  must  be  decomposed  into  first  order  differential 
equations  in  order  tc  apply  the  Runge-Kutta  method. 

Premult ipl^, all  terms  in  the  equations  of  motion 
Ea.  (3.69)  by  M  and  eliminate  the  coefficient  of  Cl.  Thus 

£  +  (al  +  bH-1£)£  +  M-11I  2.  =  M-1F ( t )  (4.1) 

where  I  is  a  unit  matrix,  and  all  matrices  are  modified  as 
discussed  in  the  sectio;.  on  static  condensation.  By  letting 

u  =  £  (4.2) 
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from  which 


a  =  £  (^.3) 

The  set  of  equations  (4.1)  is  decomposed  into  2  sets  of  first 
order  ordinary  differential  equations :s 

4  =  u  (4.4a) 

a  =  m_1F  -  M-1K  a  -  (al  +  b^_1K)u  (4.4b) 


For  the  study  of  airdropped  impact,  the  initial  displacements 

q  (usually  q  =0)  and  initial  velocity  h 

Ho  J  2°  _°  are  known,  thus 

the  initial  conditions  are: 


I  =  <10  (4.5a) 

u  =  A0  (4.5b) 

Equation  (4.4)  with  initial  conditions  as  expressed  by  Eq.  (4.5) 
can  be  solved  numerically  by  the  Runge-Kutta')  method.  Eqs. 
(4.4a)  and  (4.4b)  may  be  combined  in  a  compact  form  as: 


i  =  £(t,v) 

where  t  is  the  Independent  time  variable  and 

’a 

u 

m 

The  corresponding  initial  conditions  are: 

<1, 


(4.6) 


(4.7) 


v  = 
o 


■  o 


(4.8) 


If  the  system  has  n  degrees  of  freedom,  then  the  expression 
of  Eq.  (4.6)  in  index  notation  is: 


vi  =  Si 


(t,  v. 


Vo  ) 

2n 


1,2  - 2n  (4.9) 


The  basic  equations  in  the  Runge-Kutta  method  for  solving 
the  sat  of  ordinary  first  order  differential  equations  of 
Equations  (4.9)  are: 

v  1  (t  +  h)  =  v  i  (t)  +  £■(  ai  +  2b1  +  2ci  +  d±) 


i  =  1,2 


2n 


(4.10) 


where 

*  u  should  not  be  confused  with  a  displacement.  It  is  intro¬ 
duced  here  purely  for  convenience  and  is  as  defined  by  Eq.(4.2) 


h  =  time  step  length  of  integration 


a . 

1 

b. 
1 


c . 

1 

d. 

l 


=  h{gi  (t^  ,  ....  v2n)  ] 

=  h.[gi(t  +  ■jh,  v1  +  ^alf 


v2n  +  Ia2n) 1 


h ( t  +  ^h,  v1  + 


,  .....  v2n  +  jb2n) 1 


h[9i(t  +  h,  v.  +  ^ . v2n  +  ^2r)] 


2.  Digital  Computer  Program 

A  computer  program  code  has  been  developed  to  solve 
the  equations  of  motion  with  the  initial  conditions  of  the 
drop  impact  problem.  This  program  is  listed  in  the  appendix. 
The  required  input  data  are  nodal  point  coordinates,  mass 
distribution,  cushioning  forces,  initial  displacements  and 
velocities,  damping  coefficients,  and  member  properties  such 
as  cross  sectional  area,  flexural  rigidity  El,  torsional 
rigidity  GJ,  and  torsional  inertia  I  .  The  computer  program 
automatically  generates  the  structural  mass  matrix,  stiffness 
matrix  and  force  matrix  and  then  solves  the  equations  of  motion 
by  the  Runge-Kutta  method.  The  displacement,  velocity  and 
acceleration  at  any  coordinate  at  any  time  may  be  printed  out. 
The  printed  out  displacement  is  relative  to  the  position  where 
impact  begins.  From  these  displacements  the  relative  displace¬ 
ment  of  any  two  points  can  be  computed. 

This  program  is  primarily  developed  for  either  checking 
the  design  of  an  existing  vehicle  which  may  be  subjected  to 
airdrop  impact  or  as  a  guide  in  designing  a  vehicle  which  may 
be  destined  for  delivery  by  airdrop.  This  program  can  be  used 
not  only  for  vehicles  but  for  any  complex  mechanical  structure 
that  can  be  represented  by  a  grid  structure  model.  It  is 
expected  that  through  the  use  of  this  computational  procedure 
the  amount  of  actual  experimentation  required  for  developing 
cushioning  systems  can  be  materially  reduced.  If  elements 
of  the  vehicle  which  might  undergo  excessive  deformations 
during  impact  can  be  identified  before  the  vehicle  is  drop 
tested,  or  even  built,  the  computer  program  will  serve  as 
a  design  tool  as  well  as  an  aid  in  cushioning  design. 


FACTORS  AFFECTING  THE  RESPONSE  TO  IMPACT 

The  previous  analysis  shows  that  the  dynamic  response 
of  a  structure  is  dependent  upon  the  nature  of  the  applied 
forcing  function,  the  elastic  and  inertia  properties  of  the 
structure,  and  the  damping  characteristics  of  the  system.  By 
study  of  these  factors,  information  on  design  of  cushioning 
systems  and  the  design  of  the  structure  itself  may  be  obtained. 
Insight  into  the  appropriateness  of  the  lumped-parameter  model 
should  also  be  provided.  In  following  sections,  these  factors 
are  discussed  in  detail.  The  effects  of  the  tires  on  the 
structural  response  to  the  impact,  and  the  dynamic  behavior  of 
the  engine  on  the  rubber  supports  will  also  be  considered. 
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1.  Shape  of  Applied  Force 

In  an  airdrop  of  a  vehicle  it  is  fastened  to  a  platform 
which  holds  the  cushioning  system  in  place  but  provides  no 
support  for  the  vehic]e.  In  order  to  avoid  excessive  damage 
to  the  vehicle,  cushioning  materials  are  placed  between  it 
and  the  platform.  According  to  previous  experimental  inves¬ 
tigations,  the  force  transmitted  to  the  structure  through 
the  cushioning  material  is  essentially  independent  of  the 
degree  of  crushing.  The  amplitude  and  duration  of  the  force 
are  generally  dependent  upon  the  arrangement  of  the  cushioning 
materials  and  the  crushing  characteristics  of  these  materials. 

To  study  the  response  of  a  structure  to  the  force  applied  by 
the  cushioning,  the  response  of  the  structure  due  to  the  appli¬ 
cation  of  a  rectangular  impulse  can  be  examined.  The  effect 
of  amplitude  and  duration  of  the  force  on  the  behavior  of 
the  structure  is  cf  particular  interest.  Three  rectangular 
pulses  with  different  amplitudes  and  durations  and  one  trian¬ 
gular  pulse,  shown  in  Fig.  5.1,  were  chosen  for  this  study. 

For  simplicity,  the  area  of  each  pulse  was  kept  the  same. 

Thus  the  momentum  imparted  to  the  structure  is  the  same  for 
each  pulse.  In  order  to  simplify  the  analysis,  the  centroid 
position  of  the  pulse  has  been  chosen  as  a  characteristic 
parameter.  In  Fig.  5.1,  case  A  is  the  impulse  for  a  design 
acceleration  of  17. 5g.  Case  B  and  C  represent  the  same 
impulse  as  A,  but  with  different  time  duration.  Case  D  is 
a  triangular  pulse  with  the  same  area  and  time  duration  as 
case  a. 

Neglecting  the  damping,  the  maximum  displacements  and 
peak  accelerations  at  all  nodal  points  in  Fig.  2.1  for  all 
cases  have  been  calculated.  For  the  purpose  of  demonstration 
and  discussion,  the  maximum  displacements  and  peak  accelera¬ 
tions  at  node  18  are  shown  in  Fig.  5.2  and  5.3  respectively. 

This  node  is  on  the  frame  over  the  left  rear  wheel  of  the 
truck.  There  is  no  cushioning  force  applied  directly  at 
this  point.  Fig.  5.2  shows  that  the  maximum  vertical  displace¬ 
ment  of  this  point  as  a  result  of  the  application  of  the  four 
different  pulse  shapes  to  the  truck  is  linearly  proportional 
to  the  time  to  the  centroid  of  the  pulse  shape. 

The  peak  accelerations  produced  by  each  of  the  different 
impulses  are  shown.  These  results  suggest  that  both  the  dis¬ 
placement  and  acceleration  produced  by  a  given  impulse  depend 
essentially  on  the  time  to  the  centroid  of  the  area  under  the 
force- time  (impulse)  curve. 

2.  Structural  Properties 

The  structure  of  an  M-37  truck  is  so  complicated  that 
simplifications  must  be  made  for  analytical  studies.  If  the 
truck  is  to  be  represented  by  a  simplified  model,  the  questior 
of  how  to  estimate  the  stiffness  of  the  structure  must  be 
answered.  As  a  part  of  ihe  attack  on  this  problem  the  signi¬ 
ficance  of  changes  in  the  impact  response  with  variation  in 
the  structural  properties  should  be  investigated.  Consider 
now  the  mathematical  model  of  the  M-37  truck  shown  in  Fig.  2.1. 
Three  different  sets  of  values  of  stiffness  are  assumed: 

Case  1:  All  member  stiffnesses  are  estimated 
based  on  the  truck  frame  only 
Case  2:  All  member  stiffnesses  of  case  1  are 
multiplied  by  five 
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Fig.  5.3  Impulse  Shape  Effects  on  Peak  Acceleration 
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Case  3:  All  member  stiffnesses  of  case  1  are 
multiplied  by  ten. 

Computed  total  displacements  and  peak  accelerations 
along  the  main  frame  joint  No.  1  through  No.  10  for  the  above 
three  cases  are  shown  in  Fig.  5.4  and  5.5.  All  calculations 
are  assumed  ta  be  for  a  10  ft.  equivalent  drop  height  with 
a  design  acceleration  level  of  17. 5g  and  no  damping.  The 
centroid  of  the  honeycomb  area  is  located  intentionally  6 
inches  toward  the  front  from  the  center  of  gravity  of  the 
loaded  truck.  This  causes  a  non-uniform  crushing  of  the 
honeycomb  cushioning  as  shown  in  Fig,  5.4.  The  most  uniform 
displacement  occurs  for  case  3,  the  truck  with  the  greatest 
stiffness.  However ,  Case  2  deflection?  differ  very  little 
from  those  of  Case  3.  The  peak  acceleration  curves  shown  in 
Fig.  5.5  indicate  that  case  3  also  has  the  most  uniformity 
in  the  acceleration  at  all  points  along  the  truck  frame. 

The  differences  in  peak  accelerations  among  the  three  cases 
at  a  given  point  are  small.  All  curves  are  very  close  toge¬ 
ther.  Thus  it  appears  that  a  reasonable  impact  response  can 
be  obtained  by  using  very  rough  approximations  to  the  struc¬ 
tural  stiffness. 

3.  Damping 

The  amount  of  structural  damping  in  the  vehicle 
must  be  determined  by  experimentation.  If  experimental  data 
or  reliable  information  is  not  available,  no  damping  should 
be  assumed  in  the  analysis  since  the  omission  of  damping 
results  in  conservative  estimates  of  deflections,  and  conser¬ 
vative  cushioning  system  designs.  Numerical  computations 
for  the  M-37  model  (Fig.  2.1)  show  that  at  all  points  of  the 
truck  frame  the  absolute  magnitudes  of  displacements  are 
increasing  as  damping  decreases  but  the  peak  accelerations 
are  affected  very  little.  Fig.  5.6  shows  the  maximum  crushing 
displacements  for  a  damping  ratio  range  from  zero  up  to  0.015. 
There  is  no  change  in  the  configuration  of  the  truck  for  all 
three  cases.  The  crushing  displacement  curve  with  0.009  as 
the  damping  ratio  is  the  one  which  approximates  most  closely 
the  experimental  M-37  measurements  which  are  shown  later  in 
section  6. 

4.  The  Effect  of  Tires  on  the  Response  of  Vehicle  Body 

Since  the  vehicle  body  is  connected  to  the  wheels 
through  a  leaf  spring  arrangement,  the  magnitude  and  shape  of 
forces  transmitted  to  the  vehicle  body  as  the  result  of  an 
impact  would  be  significantly  affected  by  the  material  pro¬ 
perties  of  the  tires.  These  forces  must  be  considered  in 
computations  of  the  overall  response  of  the  vehicle  body. 

For  the  purpose  of  analysis,  the  wheel  and  tire  are 
replaced  by  a  mass,  a  spring  K,  and  da shoot  C  arrangement 
as  depicted  in  Fig.  2.3.  Using  different  sets  of  values  of 
K  and  C  for  the  tire,  the  displacement-time  curves  at  point 
No.  3  (q^)  of  the  vehicle  body,  and  the  wheel  at  point  No. 

24  (q„.)'havs  been  calculated  and  plotted  in  Figs.  5.7,  5.8, 
5.10,  and  5.11.  In  Figs.  5.r,  and  5.12,  the  maximum  displace¬ 
ment  and  acceleration  of  nodal  point  No.  3  and  the  wheel  have 
been  plotted  as  a  function  of  spring  constants  and  damping 
ratios . 

In  general,  the  larger  the  spring  constants  and  the 
damping  ratio,  the  smaller  the  displacement  of  the  vehicle. 
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Fig,  5.^  Structural  Stiffnass  Effect  on  Maximum 
Displacements. 
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Fig.  5.5  Structural  Stiffness  Effect  on  Teak  Accelerations. 


32 


<1 


fW.Vmwi  aasfcg  r.  a  i  dagAeca  i™lwfTTmJ«i 

23  4  567  8  9  10 

NODAL  NUMBER 


c  *  DAMPING  RATIO 


c  ■  0 . 015 


- - 


C  s  o  .009*' 


c= o.ooS" - 


cs  0.0 


\ 


.6  Structural  Damping  Effect  on  Maximum  Displacements, 


33 


q,  (FRAME) 


-  I 


1-  \ 


W 

Ui 

X 

o 


s 

* 

< 

_i 

ft 

o 


-2  1- 


-  3 


-4 


-6 


10 


20  30  40 

TIME  -  MILLISECONDS 


50 


\ 

\ 


C  *  0.2 


V 


\ 

\ 


\  \  ^ _ 

\ 

\ 

\ 

\ 

\ 


K  =  10,500  lb  /in 


/K  =  7000  lb /in 

/  / 

/  / 

/  / 

/  / 

/  / 

K  =  3500  lb /in 


\ 


K  =  0  lb  /  in 


Fig.  5.7  Spring  Constant  Effect  of  Tire  on  Frame 
Displacement . 
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Fig,  5.8  Spring  Constant  Effect  of  Tire  on  Wheel 
Displacement, 
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Fig.  5.11  Internal  Damping  Effect  of  Tire  on  Wheel 
Displacement . 
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Pig,  5.12  Internal  Damping  Effect  of  Tire  on  Maximum 
Displacements  and  Accelerations. 


It  is  also  seen  in  Figs.  5.9  and  5.12  that  the  wheels  have 
less  displacement  than  the  vehicle  body.  This  is  the  reason 
that  fewer  cushion  pads  are  required  under  the  tires.  The 
peak  acceleration  of  the  vehicle  body  does  not  show  any  sig¬ 
nificant  variations  as  the  spring  constant  and  the  damping 
ratio  of  the  tires  vary  in  the  range  of  this  investigation. 

5.  The  Effect  of  Engine  Supports  on  the  Behavior  of  Engine 
The  idealized  engine  supporting  arrangement  is  shown 
in  Fig.  2.2.  Computational  results  indicate  that  the  values 
of  the  spring  constant  K  and  damping  ratio  C  in  Fig.  2.2  affect 
neither  the  peak  accelerations  nor  the  maximum  displacements 
at  any  point  ;n  the  vehicle  except  the  engine  itse^r  during 
the  period  of  impact.  Figs.  5.13  through  5.16  show  the  effects 
of  the  elastic  stiffness  and  internal  damping  of  engine  supports 
on  the  dynamic  responses  of  the  engine.  As  in  the  case  of 
the  tire,  in  general,  the  larger  the  spring  constant  and  inter¬ 
nal  damping  of  the  supports,  the  smaller  the  engine  relative 
displacement  to  the  truck  frame. 

The  peak  acceleration  of  the  engine  is  not  influenced 
much  by  the  spring  constant:  of  the  supports  but  it  decreases 
with  an  increase  in  internal  dampirg  of  the  support. 

EXPERIMENTAL  INVESTIGATION  -■  M-37  TRUCK  DROP  TEST 
To  obtain  some  experimental  data  for  comparison  with 
the  computed  results  as  described  in  the  previous  chapters, 
three  drops  of  the  M-37,  3/4  ton  truck  with  a  1500  lb.  simu¬ 
lated  load  of  sandbags  have  been  made  from  a  drop  height  of 
10  ft,,  and  at  a  design  acceleration  of  17.5  g. 

The  cushioning  system  used  for  all  drops  is  shown 
in  Fig.  6.1.  Typical  design  calculations  for  such  cushion¬ 
ing  systems  are  shown  in  the  Appendix.  The  cushioning  mater¬ 
ial  used  throughout  this  series  was  80-0-1/2  paper  honeycomb 
with  a  characteristic  stress-strain  curve  as  shown  in  Fig.  6.2 
and  an  average  crushing  stress  of  6400  psf.  The  truck  was 
rigged  for  drop  by  attaching  lifting  plates  and  shackles  to 
each  of  the  wheels.  The  entire  rigging  is  she  in  Fig.  6.3. 

Accelerations  are  measured  wich  fluid  damped  resistance 
type  accelerometers.  Engine  displacements  relative  to  the 
truck  frame  are  measured  with  slide-wire  type  transducers. 

A  special  deflection  gage  shown  in  Fig.  6.4  is  used  to  measure 
the  displacements  at  other  positions  on  the  vehicle  body. 
Accelerometers  and  displacement  gages  were  mounted  on  the 
vehicle  in  the  following  positions: 

Drop  No.  1  (Series  No.  M-37-17) 

Accelerometers  at 

Engine  front  (Joint  No.  23) 

Thin  plate  above  the  winch  (Joint  No.  21) 

Rear  bumper  (Joint  No.  36) 

Right  rear  wheel  (Joint  No.  33) 

Right  middle  frame  (Joint  No.  6) 

Displacement  gages  at 

Engine  front  -  Relative  to  frame  (Joint  No.  23 
relative  to  No.  3) 

Rear  frame  -  Relative  to  wheel  (Joint  No.  18 
relative  to  No.  35) 
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Fig.  5.15  Internal  Damping  Effect  of  Engine  Support  on 
Displacements  of  Engine. 
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Pig.  5.16  Internal  Damping  Effect  of  Engine  Support  on 

Maximum  Displacement  and  Acceleration  of  Engine. 
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Fig.  6.3  Rigging  for  Lifting  the  Vehicle. 
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Pig.  6.4  Rotary  Displacement  Transducer. 
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Drop  No.  2  (Series  No.  M-37-18) 

Accelerometers  at 

Frame  (Joint  No.  2) 

Transmission  (Joint  No.  28) 

Rear  Differentia1  (Joint  No.  34) 

Displacement  gages  at 

Front  frame  (Joint  No.  3) 

Rear  frame  (Joint  No.  18) 

Transfer  case  (Joint  No.  30) 

Engine  rear  -  Relative  to  frame  (Joint  No.  38 
Relative  to  No.  4) 

Drop  No.  3  (Series  No.  M-37-19) 

Accelerometers  at 

Engine  front  and  rear  (Joints  No.  23  and  No.  38) 

Truck  bed 

Displacement  gages  at 

Front  frame  (Joint  No.  1) 

Middle  f^ame  (Joint  No.  5) 

Left  wheel  (Joint  No.  35) 

All  acceleration  and  displacement  data  were  recorded  by  a 
magnetic  tape  system  and  re-recorded  from  the  tape  to  a 
visible  record  on  paper  using  a  Visicorder  oscillograph. 

In  addition  to  acceleration  and  displacement  records, 
high-speed  motion  pictures  were  taken  for  all  drops.  These 
pictures  were  studied  to  determine  the  gross  -spects  of  vehi¬ 
cle  motion  during  impact. 

Numerical  values  used  in  the  computation  for  the  M-37 
truck  model  are  as  follows: 

Drop  height  10  ft.  (or  impact  velocity  302  in/sec.) 

Design  acceleration  level  -  17.5  g. 

Impact  duration  =  45.2  milliseconds 

Average  crushing  stress  of  paper  honeycomb  =  6400  psf. 
Dimensions  and  arrangement  of  cushioning  pads  (see  Fig. 
6.1) 

Applied  cushioning  forces: 

Concentrated  forces  - 

Wheels  and  tires  each  F  =  350  x  18.5  =  6480  lb. 
Differentials  each  F^  =  F^  =  480  x  1  8 . 5  = 

8890  lb. 

Transmission  FTr  =  200  x  18.5  =  3700  lb. 

Gear  Reducer  F  =  300  x  18.5  -  5550  11 

gr 

Uniform  distributed  forces: 

At  pad  A.,  F.  =  20600/(2  x  23)  =  448  lb/in. 

At  pad  A~ ,  F'l  =  13900/29  =  475  lb/in. 

At  pad  A^ ,  F^  =  13100/19  =  697  lb/in. 

(See  Appendix  B) 

The  weights  of  some  of  the  vehicle  components  were  obtained 
from  the  M-37  technical  manual7  and  others  were  assumed.  These 
weights  are  listed  as  follows: 

Total  truck  weight  =  5,390  lb. 

Load  weight  =  1,500  lb. 

Discrete  masses: 
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600  lb. 


Engine 

Wheels  and  tires 

(each;  350  lb. 

Transmission  200  ib. 

Gear  Reducer  3UU  ib. 

Differentials  (each) 480  lb. 

Uniformly  distributed  masses: 

All  weights  not  included  in  discrete  masses  distri¬ 
buted  as  follows: 

Longitudinal  frame  members: 

Joint  No.  1  to  No.  6,  and  No.  11  to  No.  16, 

=  570/2  x  97  =  2.94  lb/in. 

Joint  No.  6  to  No.  10,  and  No.  16  to  No.  20 
=  2560/2  x  89  =  14.38  lb/in. 

Transverse  frame  members: 

Joint  No.  2  -  No.  12  (Winch)  =  10  lb. /in. 

Joint  No.  3  -  No.  13  =2.3  lb/in. 

Joint  No.  5  -  No.  15  =1.2  xb/in . 

Joint  No.  6  -  No.  16  =0.8  lb/in. 

Joint  No.  10  -  No.  20  =  0.8  lb/in. 

Except  for  Joint  2  to  Joint  12  these  weights  were  arbi¬ 
trarily  assigned.  They  are  approximately  proportioned  to  the 
member  cross  section.  All  member  properties,  calculated  using 
truck  frame  cross-sectional  dimensions  only,  are  listed  in 
Table  6.1.  The  tire  spring  constants  are  assumed  to  be  7000  lb/ 
in  and  the  damping  ratio  0.2.  The  spring  constants  of  the 
engine  supports  are  12,000  lb/in  at  front  and  40,0C0  lb/in  at 
rear  and  the  damping  ratio  is  0.6.  The  structural  damping 
ratio  is  assumed  to  Ke  0.009.  Joint  coordinates  in  an  x-y  sys¬ 
tem  with  the  origin  at  the  rront  of  the  truck  as  shown  in  Fig.  2.1 
are  given  in  Table  6.2. 

The  measured  results  (solid  line)  are  plotted  with  the 
calculated  results  (dotted  line)  in  Figs.  6.5  to  6.20.  These 
results  are  discussed  in  the  next  chapter. 
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Start  Err1 
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I  -  in4 

X 

Torsional 

Constant 

-r  •  4  i 

J  -  in 

1 

1-6 

1.6 

4.29 

4.39 
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2.0 

10.8 

11.0 

0.027  1 

11  -  16 

1.6 

4.29 

4.39 
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0.5 

0.01 
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0.01 

0.5 

0.01 
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3.4 

4.47 

12.57 

10.6 
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40.0 

30.0 

60.0 

60.0 

33  -  35 

40.0 

30.0 

60.0 

60.0 

4  -  14 

100.0 

100.0 

100.0 

100.0 

27  -  28 

100.0 

100.0 

100.0 

100.0 

5  -  15 
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1.4 

0.017 

6  -  16 

1.3 
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1.4 

0.017 

7-17 

2.0 
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8.0 

14.4  j 

9  -  19 
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2.5 

5.2 

0.024 

10  -  20 

1.8 
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5.2 

0.024 

3-4 

13  -  14 

7-9 

2.0 

0.167 

0.92 

0  917 

l7  -  19 

_ 

Table  6.1  Member  Properties 


Joint 

No. 

Coordinates 

Joint 

No. 

Coordinates 

X  -in. 

Y  -in. 

X  -in. 

Y  -in. 

1 

0 

15.0 

21 

6.0 

0 

2 

6.0 

15.0 

22 

17.0 

0 

3 

17.0 

15.0 

23 

17.0 

0 

4 

51.0 

15.0 

24 

34.0 

15.0 

5 

87.0 

18.6 

25 

34.0 

5.0 

6 

97.0 

19.0 

26 

34.0 

-15.0 

7 

117.5 

19.0 

27 

5x .  0 

0 

8 

147.5 

19.0 

28 

66.0 

0 

9 

177 . 5 

19.0 

29 

87.0 

6.0 

_0 

186.0 

19.0 

30 

87.0 

-6.0 

11 

0 

-15.0 

31 

97.0 

6.0 

i2 

6.0 

-15.0 

32 

97.0 

-6.0 

13 

17.0 

-15.0 

33 

147.5 

19.0 

14 

51.0 

-15.0 

34 

147.5 

5.0 

15 

87.0 

-18.6 

35 

147.5 

-19.0 

16 

97.0 

-19.0 

36 

186.0 

12.0 

17 

117.5 

-19.0 

37 

186.0 

-12.0 

18 

147  5 

-19.0 

38 

51..  0 

15.0 

19 

177.5 

-19.0 

39 

51.0 

-15.0 

20 

186.0 

-19.0 

Table  6.2  Joint  Coordinates 
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6.5  Measured  and  Computed  Acceleration  of  q_  -  winch 
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Pig.  6«fi  Measured  and  Computed  Acceleration  of  qP,  -  winch 
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Fig.  6.7  Measured  and  Computed  Acceleration  of  q,  -  Middle  Frame 
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6.9  Measured  and  Computed  Acceleration  of  q?R  -  Transmission. 


COMPUTATION 


59 


ransfer  Case 


COMPUTATION 


62 


Fig.  6.15  Measured  and  Computed  Relative  Displacement  of  (q„_  -  q-,)  -  Engine 
Relative  to  the  Frame.  ^  ^ 
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Plg‘  6’16  Relative  °f  '^8  '  Q*>  -  Engine 


17  Measured  and  Computed  Displacement  of  -  Front  Frame 


Fig.  6.19  Measured  and  Computed  Displacement  of  R  -  Rear  Frame. 
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DISCUSSION 

L.  Mathematical  Model 

The  analysis  oi  the  complex  structure  presented  in 
this  work  is  based  on  the  possibility  of  modeling  the  struc¬ 
ture  by  a  lumped-parameter  system.  Modeling  of  the  structure 
is  somewhat  arbitrary.  In  general,  however,  the  following 
two  faccors  need  to  be  carefully  considered. 

(a)  Select  and  sort  out  the  components  in  the  struc¬ 
ture  whose  motion  is  to  be  studied.  It  is  impractical  to  try 
to  include  all  the  components  of  the  structure  in  a  model. 

It  is  necessary  for  the  designer  to  select  the  components 
whose  motions  are  considered  to  be  most  important.  For  example, 
in  an  airdrop  of  a  vehicle,  the  motion  of  the  engine  block, 
transmission,  differentials,  and  wheels  may  be  the  important 
parameters  for  the  proper  design  of  a  cushioning  system  to 
protect  the  vehicle.  Consequently  the  model  should  be  so 
designed  as  to  represent  the  motions  of  these  components  as 
realistically  as  possible.  On  the  other  hand,  if  the  task  is 
to  design  a  cushioning  system  for  the  protection  of  a  car 
radio,  a  completely  different  model  will  be  needed.  The  mode] 
shown  in  Fig.  2.1  is  developed  for  the  purpose  of  designing 
a  cushioning  system  for  the  protection  of  the  vehicle  during 
airdrop.  The  weights  of  the  engine,  transmission,  transfer 
care,  differentials,  and  wheels  are  modeled  by  concentrated 
masses  at  proper  nodal  points.  The  winch  and  the  frame  of 
the  truck  were  replaced  by  two  longitudinal  beams  coupled 
together  by  eight  transverse  beams.  The  weight  of  the  truck 
frame  and  the  load  are  further  assumed  to  be  uniformly  dis¬ 
tributed  along  the  two  longitudinal  beams.  The  stiffnesses 
of  all  the  connecting  members  in  Fig.  2.1  are  calculated  on 
a  static  basis  using  the  dimensions  of  the  truck  members.  It 
should  be  noted  that  Fig.  2.1  is  a  particular  model,  many 
other  models  may  be  developed  for  the  same  purpose.  The  suita¬ 
bility  of  the  present  model  will  be  discussed  in  the  section 
where  the  comparison  is  made  between  the  measured  and  computed 
results . 

(b)  Estimate  the  motion  of  the  structure.  In  model¬ 
ing  a  complex  structure,  information  on  the  environment  that 
the  structure  is  likely  to  be  subjected  to  is  also  very  impor¬ 
tant.  In  airdrop  of  a  vehicle,  the  whole  structure  is  expected 
to  land  with  an  impact  velocity  of  approximately  30  f ps ,  and 
the  rise  time  of  the  impulse  imparted  to  the  vehicle  Is  of  the 
order  of  milliseconds.  These  are  important  data.  A  continuous 
elastic  structure,  such  as  a  vehicle,  has  in  principle  an 
infinite  number  of  modes  of  vibration  that  can  be  excited.  A 
lumped  parameter  system  can  be  excited  in  only  as  many  modes 

as  it  has  degrees  of  freedom.  Therefore  a  lumped  parameter 
system  can  not  be  expected  to  accurately  represent  a  continuous 
structure  under  conditions  in  which  the  higher  modes  of  the 
structure  might  be  excited.  The  long  rise  cime  of  the  impulse 
applied  to  the  vehicle  means  that  the  higher  modes  will  not 
be  excited  and  the  lumped  parameter  model  c  n  therefore  pro¬ 
vide  a  suitable  approximation  to  the  motion  of  the  prototype 
system.  On  the  other  hand,  an  input  force  with  a  very  short 
rise  time  would  tend  to  excite  higher  modes  of  vibration.  As 
a  consequence  a  lumped  parameter  model  would  not  be  able  to 
represent  the  motion  of  the  prototype  structure. 
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In  the  present  work,  the  whole  truck  body  is  modeled 
by  a  lumped  parameter  system  with  39  nodal  points.  Since 
there  are  six  degrees  of  freedom  at  each  nodal  point,  the 
total  number  of  degrees  of  freedom  of  the  system  is  234. 

Only  a  normal  impact  (in  the  vertical  direction)  is  considered 
here,  hence  the  vertical  displacements  of  the  nodal  points 
are  of  the  main  concern.  The  six  displacements  possible  at 
each  nodal  point  consist  of  three  translational  motions  and 
ihree  rotational  motions.  Two  of  the  translational  displace- 
*nents  and  one  of  the  rotational  displacements  can  be  disre¬ 
garded.  Thus  the  total  number  of  degrees  of  freedom  are 
reduced  to  117.  This  number  can  be  further  reduced  to  39 
by  eliminating  the  remaining  two  rotational  components  at 
the  nodal  points  by  the  static  condensation  process.  4 

In  the  numerical  computations,  the  element  mass  matrix 
and  stiffness  matrix  are  calculated  according  to  equations 
(3.19)  and  (3.29)  respectively.  The  element  end  forces  are 
calculated  by  equations  (3.31)  and  (3.32).  Only  concentrated 
forces  and  uniformly  distributed  forces  are  considered  in  the 
analysis.  The  direct  stiffness  method  is  used  to  formulate 
the  system  mass  matrix  and  stiffness  matrix.  The  structural 
damping  matrix  is  treated  as  a  linear  combination  of  the  mass 
and  stiffness  matrices  of  the  system.  The  amount  of  computer 
storage  required  by  the  program  for  this  analysis  is  30,000 
words.  For  the  CDC  660  j  computer,  the  formulation  of  the 
equations  of  motion  requires  three  minutes.  The  numerical 
integrations  by  the  Runge-Kutta  method  require  1.5  minutes 
for  an  impact  duration  of  45  milliseconds,  using  an  inte¬ 
gration  step  size  of  0.1  milliseconds. 

2.  Comparison  of  Measured  and  Computed  Results 

From  Fig.  6.5  to  6.20,  it  may  be  seen  that  agreement 
between  measured  and  computed  results  is,  in  general,  not 
very  good,  so  far  as  the  shapes  of  the  curves  are  concerned. 
However,  the  agreement  in  the  amplitudes  of  displacements  and 
accelerations,  with  the  exception  of  the  accelerations  of 
wheels  and  differential  is  quite  good.  The  maximum  discre¬ 
pancy  between  measured  and  computed  displacements  occurred 
at  che  front  iram  and  is  18%.  The  maximum  discrepancy 
between  measured  and  compui.  accelerations  at  the  engine 
is  34%.  For  the  acceleration  ot  the  differential  the  maxi¬ 
mum  discrepancy  is  70%,  and  for  the  wneels  the  maximum  dis¬ 
crepancy  is  50%.  The  displacements  are  believed  to  be  more 
significant  as  indicators  of  possible  damage,  than  the  acceler¬ 
ations.  Hence  it  is  encouraging  to  find,  in  this  example, 
that  the  displacements  can  be  computed  with  the  degree  of 
success  indicated.  The  actual  shapes  of  the  displacement  and 
acceleration  curves  are  probably  of  little  consequence  except 
for  the  indications  they  give  of  the  adequacy  of  the  lumped- 
parameter  model. 

The  agreei  between  computed  and  measured  displace¬ 
ments  and  accelerations  may  be  improved  in  several  ways, 
namely : 
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(a)  Change  uhe  centroid  position  of  the  impact 

forcing  function  by  adjustments  in  the  cusi  jn- 
ing  areas. 

( b ’<  Adjust  the  structural  damping  ratio. 

(c)  Adjust  the  spring  constants  and  damping  ratios 
for  a  more  accurate  representation  of  rubber 
properties . 

(d)  Replace  the  constant  cushioning  force  applied 
to  the  wheel  by  a  linearly  varying  force. 

Finally,  it  should  be  noted  that  computed  and  measured 
results  should  not  be  expected  to  agree  in  all  details  since 
the  model  on  which  computations  are  based  is  only  an  approxi¬ 
mation  of  the  real  structure,  and  the  initial  conditions  used 
for  computations  are  idealized  approximations  to  the  initial 
conditions  that  actually  occurred  in  the  experiments.  In 
general,  in  lumped-parameter  modeling,  the  parameters  can  be 
adjusted  to  provide  agreement  between  computed  frequencies 
and  corresponding  frequencies  in  the  prototype  structure,  or 
the  adjustment  can  be  made  so  amplitudes  of  displacements  and 
accelerations  at  selected  points  agree  between  model  and  pro¬ 
totype.  However,  agreement  between  frequencies,  and  between 
amplitudes  cannot  both  be  achieved  simultaneously.  A  decision 
must  be  made  before  parameter  values  are  assigned,  as  to  what 
quantities  are  to  be  matched. 

CONCLUSIONS 

From  the  previous  discussions  the  following  conclusions 
are  believed  to  be  justified. 

a.  Displacements,  velocities  and  accelerations  at 
various  points  of  a  complex  mechanical  structure  can  be  satis¬ 
factorily  predicted  using  a  lumped  parameter  mathematical 
model  and  a  numerical  computation  procedure. 

b.  The  dynamic  response  of  a  complex  structure  when 
properly  cushioned  and  subjected  to  an  impact  loading  is  not 
sensitive  to  the  elastic  properties  of  the  structure.  Thus 
the  elastic  coupling  between  masses  in  the  lumped  parameter 
model  need  not  be  known  to  any  great  degree  of  precision. 

c.  Structural  damping  dissipates  a  considerable 
amount  of  energy,  and  as  a  consequence  decreases  the  displace¬ 
ments.  However,  the  peak  accelerations  a:  various  points  in 
the  structure  are  affected  very  little  by  the  structural 
damping . 

d.  The  forcing  function  is  the  major  factor  which 
affects  the  dynamic  response  of  the  system.  It  must  be  repre¬ 
sented  as  exactly  as  possible.  This  means  that  the  cushioning 
characteristics  are  very  important. 

e.  The  procedure  for  handling  the  engine  and  transfer 
case  in  the  analysis  can  be  applied  to  any  rigid  discrete  mass 
which  cannot  be  included  in  the  elastic  properties  of  the  struc¬ 
ture  . 

f.  The  developed  computer  program  can  be  employed  to 
predict  the  dynamic  response  of  any  complex  mechanical  struc¬ 
ture  if  the  structure  can  be  represented  by  a  grid  type  model. 
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This  program  can  be  used  to  determine  where  large  relative 
deflections  may  be  expected  to  occur  if  an  existing  struc¬ 
ture  is  subjected  to  an  airdrop  impact,  or  to  pinpoint  trou¬ 
ble  points,  so  far  as  airdrop  is  concerned,  in  vehicles 
still  in  the  design  stage. 

RECOMMENDATIONS  FOR  FURTHER  STUDIES 

a.  Only  vertical  motion  was  considered  in  this  study 
This  is  an  oversimplification  of  the  real  situation  because 
non-vertical  and  con-planar  motion  is  always  possible  due  to 
wind  drift  and  system  oscillation  in  the  airdrop  process,. 

The  effects  on  the  dynamic  response  a  structure  due  to 
non-vertical  and  non-planar  impact  juld  be  investigated. 
This  can  be  done  by  introducing  six  agrees  of  freedom  at 
each  nodal  point  of  the  lumped  parameter  model  in  the  mathe¬ 
matical  analysis.  However,  shearing  properties  of  cushioning 
materials  will  be  needed. 

b.  For  the  purpose  of  structural  design,  the  compute 
cion  of  stresses  in  all  members  should  be  included  in  the  com 
puter  program.  It  may  be  necessary  to  consider  the  behavior 
of  elast-i-plastic  structural  elements  when  subjected  to  high 
velocity  impact.  However,  the  step-by-step  numerical  integra 
tion  method  used  in  this  study  can  be  modified  to  include  the 
yield  conditions  of  the  element. 
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LIST  OF  SYMBOLS 


(t) 


■+  f 

,j.y- 


(x  t) 
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Transformation  matrix  relating  W  and  U. 

Matrix  relating  stress  and  strain  S  =  A  C 
Unit  vectors 

Matrix  relating  strain  and  displacement  e  =  b  U 

Transformation  matrix  c[  =  B  U 

Damping  matrix  in  u  coordinates 

Damping  matrix  in  q  coordinates 

A  damping  coefficient 

Constants 

Damping  force  matrix  D  =  -C  2. 

Position  vectors 
Strain  matrix 
Bending  strain 
Torsional  strain 
Force 

Cushion  force  matrix 

Equivalent  concentrated  forces  acting  at  element  ends 
due  to  distributed  forces. 

Force  vector  in  global  coordinates 

Shear  Modulus 

Torsional  inertia 

Rotatory  inertia 

Unit  vectors  along  datum  coordinates 
Torsional  constant  (polar  moment  of  inertia) 

Element  stiffness  matrix  in  u  coordinates 

Element  stiffness  matrix  transformed  to  global  coordinates 
Stiffness  matrix 
A  spring  constant 
Element  length 

Element  m.^s  matrix  in  u  coordinates 

Element  mass  matrix  transformed  to  global  coordinates 

Mass  matrix 

Moment 

External  force  matrix,  acting  as  an  element 
Joint  coordinates 
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Displacements  in  direction  of  datum  coordinates  X,Y,Z 
r  =  X/L  Nondimensionai  parameter 
s  =  Y/L  Nondimensionai  parameter 
S_  Stress  matrix 

Bending  moment 
S  Twisting  moment 

t  =  Z/L  Nondimensionai  parameter 
T  Kinetic  energy 

U  Element  displacements  referred  to  local  coordinates 

x,y ,  z 

U  Displacement  matrix 

V  Volume,  strain  energy 

W.  Displacements  referred  to  element  coordinates  x,y,z 

W  Work 
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APPENDIX  A 


low  Chart  and  Computer  Program 


C  NUMERICAL  SOLUTION  OP  A  LUMPED-PARAMETER  MODEL  FOR 
C  AIRDROPPED  STRUCTURES  b*  USING  RUNGE-KUTTA  METHOD. 

C  CDC  6600  ,  FORTRAN  IV,  STORAGE  REQUIRED  30000  WORDS 

PROGRAM  ARDROP  ( INPUT , OUTPUT ) 

COMMON  N,$(120,120),SMK(40,40) 

COMMON  TM(20)  ,CP(20),SF(1?0),WG(120)AN1,LN2,SLOP 
COMMON  NDP,JP1 (50) , JP2C50, ,DP(50) 

COMMON  GAMA, SK, BETA  , MTF 

DIMENSION  W(6,6),RC6,6),WM(6),WF(6),FD(6),EM(6),AB(6), 

1  DA(6,6),DAT(6,6),tD(6,6),A<6,6),fl(6,6),D(6,6), 

2  SK(120)»SM(120)»X(43),Y(40),MR(120),TR(120),NR(120), 

3  MV(120),MO(120),NQ(1?0),YJ(80),DF (5,80)  ,SF0(120) 

4  ,  EF (6) , SD( 120) 

PRINT  6001 

C  READ  IN  JNT  »N0.  OF  JOINTS 
C  MT  oNO.  OF  ELEMENTS 

C  LAGMsNO.  OF  DISCRETE  MASSES  SUPPORTED  BY 

C  POINTS  MORE  THAN  TWO 

C  IQR  B-l,  ELIMINATING  ALL  ROTATIONAL  COORDINATES 

C  *  0,  NO  REDUCING  IN  COORDINATES 

C  «N0.  OF  COORDINATES  WILL  BE  ELIMINATED 

READ  5040,  JNT , MT , L AGM , I  OR 
PRINT  5040,  JNT,MT,LAGM, IQR 
C  --  READ  IN  EsMOOULUS  OF  ELASTICITY 
C  GbSHEAR  MODULUS 

READ  5010, E,G 
PRINT  5010, E,G 
IR=3*JNT 

00  10  1  =  1, IR  S  SF(I)«0.  5  W  G ( I ) =  0 .  $  SK ( I ) *0 , 

SO ( I ) =0 . 0  %  SM(I)30, 

DO  10  J=1,IR  5  S(I,J)B0, 

10  CONTINUE 

00  20  1=1, JNT 

C  ---  READ  In  J  a  JOINT  NUMBER 
C  X(J) , Y(J)sX,Y  COORD.  OF  NODE  J 

C  WM(1),WM(2),WM(3)bMASS,MA3S  MOMENT  OF  INERTIA 

C  ABOUT  X  AND  Y  RESPECTIVELY  AT  MODE  J 

C  WF ( I ) aFORCES  APPLIED  AT  NODE  J 

READ  5020,  J,X(J),Y<J:,(wm(K),K«1,3),(WF(K),Kb1,3) 

PRINT  5020, J,X(J)  ,Y(J, , (WM(K),K*1,3)»(WF(K),K*1,3) 

K 1= J*3-2  $  WG(Xl)sWG(Kl )-WM(t )*386,4 

DO  15  K  =  1 , 3  S  SF  <K  1 )  wSFCKl  )  4-WF  (K)  3  8 M  ( K 1  )  =  S M  ( K 1  )  +  W K  (  K ) 

K 1 *K 1 + 1 

15  CONTINUE 

20  CONTINUE 

DO  220  IMB«1,MT 

C  ---  READ  IN  N J 1 , N J  2s JO  I  NT  NUMBER  AT  ELEMENT  ENDS 
C  MP(I)*1,  RELEA3E  THE  MEMBER  CONSTRAINT  IN  THE 

C  MEMBER  COORD .  U ( I ) 

C  =0, NO  MEMBER  CONSTRAINT  IS  RFlEASED  IN  U(I) 
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c 

c 

c 

c 

c 

c 


&  3 


26 

30 

C 

11 


<10 


60 


KF  «('  ,  NO  MEMBER  FORCES 

*i,  uniform  member  forces 

■2f  concentrated  force  at  the  member 

RO*  MAHS  PER  UNIT  LENGTH  OF  THE  MEMBER 
AA  /  A  I,  yj.G.Ts  ARS  A,  MOMENT  OF  INERTIA  ,  POLAR  MOirNT 
OF  INERTIA  AND  TORSIONAL  CONSTANT  OF 
THE  MEMBER  RESPECTIVELY 

FD  ft  UNIFORM  FORCE  FOCI)  FROM  FDC2)  TO  FD(3>,i)R 


CONCENTRATED  FORCE  AT  FD(2) 

",J1  f M»'TI . Tsj ,65  f Kp, aj,gj,RO, (FD(I) ,!•! ,3) 


READ  5030,  NJ1,NJ2, , IB- 
READ  5010,XJ,AA 
PRINT  5030,  NJ1,NJ2, (MR(I),Iul, 6), KF,AI,GJ,RO, (FOCI), 1*1,3) 

°RINT  5070, XJ, AA 

Xi  =  X(NJ2)-XfNJn  $ 

3  L  a  5  0  »  T  (  X  1  X  1  *  ?  1  *  Y  H 
I  r  ;  SL  5  25,25,33 

K1=NJ 1*3*2  S  «2*NJ2*3-2  $  SK(K1  )  *SK(KJ  )  +  AI 

SKCK2  3  *SK(K2)  +AI 
S ( K 1 , K 1 ) *SK  ( K  1 3  S 
IF(K1,GT,K2)  C-0  TO  26 
SCK1,K2)*S(K1,K2)*AI  5  GO  TO  220 
S(K2,K1)*S(X2,K1)»  AI  $  GO  TO  220 

Glsxi/SL  5  G2SY1/8L  S  C9«E*AI/SL**3 
***  TRANSFORMATION  MATRIX  (DIRECTION  COSINES) 

DO  11  1*1,6  $  DO  11  Ja 1 , 6 


Y1bY(NJ2)-Y(NJ1) 

SK  CK  1  ) 
S(K2,K2)*8K(K2) 

GO 


DO  11  1*1,6 
OA(I, JJB0.0 
DA ( 1 , 1) *1 ,0 
0 A ( 3 , 2 ) **G2 
D  A ( 5 , 5 ) *G 1 
DA (6,6) SGI 
***  ELEMENT 
N  (  i  ,  1  )  e  1  2 . 0 
H(fl,  1 ) a* 1 2 , 


*** 


S  DA (2, 2) *G1 
$  DA (3, 3) *G1 
S  D  A ( 5 ,6 ) *  G  2 


S 


-"KR,  Ua-12, 
W(2,2)*G*GJ/(SL*CS) 

w(a,2)*a, 
W(3,3)*«.*SL**2 
w(6,3)*2,*SL**2 
W  (b,  <1)  *6.*SL 
rf(6,6)s<l,*8L**2 
DO  <10  J*  1 , 6  $ 

WCJ,I)BW(I, J) 
IF(RO.LE,0,) 
CMsRO*SL/«20. 

R ( I , J ) *0 , 
WMC13*PO*366.<l 
CALL  EMOFOC  (EM, 
A0*AI/AA  $  A1*A0/SL 
S2*SL*SL 


STIFFNESS  *** 

5  W(2,l)«0. 

*  H (5 , 1 ) *0 , 

S  W(3,2)*0. 

S  H(5,2)*-W(2,2) 
S  >I(«,3)b6„*8L 
WU,4)«12,  S 
S  W(5,5)«W(2,2) 


S  D A (2, 3) *G2 
S  DAn»,«)*l,0 
S  DA(6,5)S-G2 

$  W(3,1)*-6,0*SL 
S  W(6, l)a-6,0*SL 

$  W (6, 2) *0, 

S  H(5r3)«0. 
W(5,fl)»0, 

$  W(6,5)*0, 


DO  <10  I* J, 6  S  W(I, J)*CS*W(I, J) 


GO 

3 


TO 

DO 


70 

60 


1*1,6  S  DO  60  J*l,6 


S  WM(3)*SL 


S  WM  C 2 ) *0 , 

1 r  WM, SL) 

J  A2*A l/SL  S  AJaXJ/AA 


78 


BEST  AVAILABLE  COPY 


i&ti  r* 


Htib  Wy  : 


75 

70 


82 


92 

110 

120 

121 

1 10 

111 

160 

170 


190 

208 

205 


2  23 


2«SS 

253 


* 


.vsSs^ijSi-i <sti 


*  *  ♦  ?■  L  f  u  E  N  T  ••  «  «  3  **  i  f  6  1 1 

*  •  • 

'  i  , :  1  *  1  c  *■  ,  *  5  F  *  .  *  *  2  % 

«(1  ,3)*-22 

“><»<*'*  '-,•50'*, *A?  I 

*M  ,61*1  3. 

•*(2,2;«160.*4J  j  9(2, 81*70.  .*J 

7  .  3  ■  »  ti ,  *  S  2  ♦  ?  6  ,  *  A  0  | 

*(3,«)»-l 3 

*(?,(■'*«  7  ,*$2-13, *A0 

w(fc,<»lr556#*5f<lt*A2  1 

• (A»6)*22, 

* ( 5 , 5  7  * i LP  .  ** J  J  B(6, 

61*4, *32*56 

1 


■  6 


S  DO  75  J»l,t> 


H  »  0 


*c:#ji*p<:,j)*c« 

»  <  J  ,  n  *  «?  r  T .  J ) 

C'On  t  I  nuE 

IF  (*.*  ,L?  »C*-3  GO  70  0P 

C*Ll  EnOFCC (EF,<F ,FO,  SL) 

DO  160  1*1,6 

I*  <*<?(!  1.LE.0)  GO  TO  160  t 
IF(w(I,IJ,CCt0»)  CO  TO  160 
OD  110  J«1,0  S  "**6l  S  IF ( J.EQ.  I)  60  TO  6* 

TB(  m)  s-w  (  I,  J) /k  (  J  ,  I  )  $  N»(M)*J  I  GO  TO  U9 

«*M-  J 

CONTI  5>J* 

CALL  C0NMS6  (*,*>,*(1,5,  I,T«) 

IECKF)  138,150,120 
00  121  6*1, M  t  K 1 »NB ( k ) 

£  F  (F1)*£F(K1 )  ♦  £  F (I)*TR(K) 

EF(ij«e, 

CONTINUE  *  IF  (RO,LE,0)  60  TO  160 

CALL  C0n*S6  (9,6, NR, 5, I,TR) 

00  131  **1,M  S  *  1 *NR ( K) 

EH(Kn«EM(Kl)*E«(n*Tft(K) 

EH(I).0, 

CONTINUE 

<l*NlJl»3-3  S  «2*NJ2#3-5  S  00  170  I»M  t  Jal*3 

**»(I)«Kl6l  t  *R(J)aK26l 

CONTINUE  I  CALL  TRAN3F(DA,DAT,0,6) 

CALL  **UL  t**DA,ED,6,6,6)  1  CALL  *»«UL  (0  AT  ,  EO .  W ,  6, 6 , 0) 

CALL  AOD*S (M, 6, S, SK, IR,HR, 1 ) 

IF  (FF  ,LE.B)  GO  TO  200 

CALL  *Ull  CDAT,EF, AB,6)  IS  00  190  1*1,6  S  K*NR(I) 

SF(p<)*8F  (k.)6AB(I) 

IFtRO.LE.0)  GO  TO  220 
CALL  *UL1 (OAT,€H, A8,6) 

00  205  1*1,6  S  F*MR(I)  f  *6 (K ) aNS (K ) -AB ( I ) 

CONTINUE 

Call  m*UL  (*#  DA, ED, 6,6,6)  S  call  M*Ul (DAT, EO, R, 6, 6 , 6) 

CALL  AOONS(R,6,S,3H, IR,HR,?) 

CONTINUE 

IJ«0  S  IKs0  S  00  290  I* 1 , I R  S  IKaIR*l 
IF (  5K ( I) ,GT,0.  J  60  TO  240 

IK«IK-1  j  IJbIJM  S  H  V  ( I J )  a  I  S  60  TO  2Si 

NQCI*)*I 

CONTINUE 
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BEST  AVAILABLE  COPY*! 


i-jittpm 


28  3 


i  ■ 


" k  <  ■  1  •  t  j 
R  F  M  0  V.  8  (  5  ,  s  h  ,  I  p  ,  m  V  ,  I J  ,  2) 

SfcMOr  C  S  f  :*  !  3  ,  u  V  ,  I  J  )  s  CALL  R  E  M  D  ►  t  WG  ,  I R ,  HV ,  I J  ) 
-  1  J 


°i  t-  . 


‘ . L  F 

I  *  i  9  ..  w  o  ' 

=  ?. 

L  ’*  :  '  '  • ' 


283  J»l,3  t  00  280  K*|,3 


:  irs  must 
M  °  r ~  k !  r  ")  E  *  j  v  3  £  c, 
*C.  .•  r  :•*.  CENTROID  OF 
<  A  •  '  '  "  i  nt ■.  >  S  In  v 

AYES  OF  T-lE 

WF ( 1 ) ,  wF<?) , Wp(3)« 


■'-OSTIS  OF  4  DISCRETE  MASS 
8E  -opATER  T*AN  TWO! 

•) r  5'jton9T  j 

“&S8 

3  E  £  5  BFTMEEn  T*E  PRINCIPAL 
.  S  S  * k-  0  X  -  a  *  1 3 

flPLlEO  MOMENTS  ABOUT  X,Y  AXES 
INDIRECTION  RESPECTIVELY 

INERTIA  ABOUT 


Rli  C 

PRINT 
READ 
PRINT 

X AaXA*3 
IFCLMS. 


-'.NO  FORCE  IN 

R C 1 , i ; , R C 2,25 , R (3, 3 )»MASS  MOMENTS  OF 
v.v  AXES  A  NO  RESPECTIVELY 

;OA10,L'JS.-  'MR(J)  ,  J*l,lMo| 

5090, L^S, (HR(JJ,J»i#LhS) 

50«2,XC,YC,XAfYA,(WFCJ),J*l,3)»(R{J,J)#J«l#S) 
50fl2,XC»YC,XA,YA,<WF(jj,J*l,3)»(R(J,J)#J«l,3) 
,1918/180,  S  YA*YA*3.IU1<»/1«0. 

LE , 3 )  00  70  310  . . 


M  N  t  L  N  S ■ 3  S  X  l  »HR  ( l )  s  K2«MR(2)  S  K3«Mpe3) 

AlrXC<2)-X(Kn  S  A2*YCK2)-Y(K1)  I  B 1 «X ( K3 ) -X (K I) 
B2*Y(K3)-Y(Kn  S  •  CC*A i *B2-A2*B  1 
00  285  J=l,3  S  KIsMRCJJ*3-2 
285  NRC J)sMACH(NO, IR,KI3 

oo  3ea  j»a,Lk'S  $  ji*j-3  s  kk*mr(jj  s  ci«x(kk)*x(kd 
C2*v(Kk)«Y(K1)  S  AA*Bl*C2-B2*Cl  S  3B« A 1 *C2-A2*C 1 
TR(11*(AA-BB*CC;/CC  $  TR(2)s-AA/CC  8  TR(3)«BB/CC 

M  :<  5  <  x  *  3  -  2 

uV(Jl)cMACH(N3,IR,KK)  J  LMK»HV(J1) 

CALL  CONNS (S,SK, I R , NR , 3 , LNK , TR , 1 ) 

CALL  C0NMS;S,5M,  IR,NIR,3,LNK,TR,2) 

CALL  C0NF(3F,IR,NR,3,LMK,TR)  S  CALL  CONF ( WC , I R , NR , 3 , LMK , TR ) 

300  CONTINUE 

CALL  REMOMS CS,SK, IR,MV,MN,  1  ) 

CALL  RENONS  CS,3m,  IR,MV,'<N,2) 

‘■■LI.  REM0F(3F-,IR,NV,m>jj  s  CAlL  REMOF  (  WG ,  I R ,  MV  ,  MN  ) 

CALL  SE3NCN3,  IR,mv,MN} 

310  CONTINUE 

A(l# IJbCOS(XA)  S  A(1 ,2)«SIN(XA)  S  A  ( 2  ,  l  ) *005 ( Y A ) 
A(2,2)bSIN(YA)  S  DO  320  J*l,3  S  K«MR(J5 


0  ( J 

, 1 >*X(K) 

-XC 

320 

D(J 

,2JeY(«T 

-YC 

00 

330  J*i, 

2  S  CO  333 

K»  1 , 3 

330 

8  (* 

,J)*A(J, 

n*DCK,2W(J, 

,2)*D(K, 1) 

DO 

390  K*l, 

3 

390 

3  (K 

,  3 )  « l , 

So 


CAL.  MlV(B,EO,M  J  CALL  ‘■<yJl('R,F0,4,l,3, 
CALc  TRANS?!  ED  *  D  *  3  *  3  3  S  CALL  MU.,L(D,8,R 
CALL  mjL1  ID,  *F  ,  AB,  3) 

EF(ns»(  3. 33  *386. 4  S  EF ( ? ) *? ,  «  FF(33*0 

CALL  u  LI (r.EF.FC.i) 

DO  370  ! 131,3  S  DO  350  J*1,IR 
* I=yR ( I  1 3 *3-2  f  I T  f  NQ(  J)  .FQ.KI  )  GO  TO  360 
350  CONTINUE  S  GO  TO  370 
360  MV ( II 3=J  S  S  F (J3»SF(J3*Ao(II) 

*G(  J3**G(J)-»FD(  II 3 
370  CONTINUE 

CALL  APOMS(R,3,S,Sm,IR,mv,23 
380  CONTINUE 
385  CONTINUE 

IF(IQR)  400,408,388 
388  CON- T  I  NUE 

C  —  IF  IQR  GREATER  THAN  ZERO 
C  READ  IN  MQ(l)sC0ORD.  TO  BE  ELIMINATED 

READ  5040,  (MQ(I) , 1*1 , IQR)  S  IQ*1 

PRINT  5040,  (MQ(I),I«1,IQR) 

GO  TO  406 


400 

continue 

I QP* JnT  *2  S  DO 

405  1*1, JNT  S  12*3*1 

K  2c  2  *  I  S  K 1 s  K  2- 

1  S  MQCK1)*I1 

405 

MG(K23*I2 

408 

I Qc 1  %  J 1 3  0 

00  60?  1*1, IP  r 

11*11+1  f  K*NQ ( 1 5 

410 

IF(IQ.GT.IQR)  GO 

TO  450 

IF(MQCIQ) ,NE.<  5 

GO  TO  430 

IQcia+1  5  GO  TO 

470 

410 

IF(40CIQ5  ,GT,k) 

GO  TO  450 

5  GO  TO 

410 

450 

IF(SM(I),GT,0.) 

GO  TO  600 

470 

JlsJl  +  1  S  M  V  ( J 1 ) s  I 

510 

MB0 

OG  560  y  1  *  1 , 1  R 
IF(SK(M1) .EQ.0.5 

M»M*  1 

GO  TO  560 

IF(mi-I)  520,515 

,540 

515 

mbm-i  %  GO  TO  560 

520 

3a«5(Hi,I5  S  GO 

TO  550 

540 

88*3(1, Ml) 

550 

TR  C  M)  a«  tjri/SK  ( 1 5 

560 

continue 

CALL  CONMS(S,SK,I 

f.  ,  MR  *  M  ,  I  ,  T R  ,  13 

CALL  CONFtSF, IR,MR,M,  I  ,TH) 

XF(SM(I)#EQ.0.)  GO  TO  6~? 

CALL  CONMS(S,SM,:.R,MR,Mf  I  ,  TR  ,  2 ) 
CALL  CONF(*G, IR,mP,m,  I  , TR 5 
600  CONTINUE 
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«9?  'j  ,*w  a  $  ,  <•*?<(!) 

SIn«J®  J  AvG**3UP/3I^ 

*0$  C  C  »-  T  T  *■<  „■  Z 

: - »C.O  IN  *C»9E*  l3*C!n;  CtS!  NUMBER 

•■:><**  >0.  Of  VISCOUS  DAMPINGS 

3  AMA*ST»UCTuRM,  DAMPING  RATIO  OP  CRITICAL  VALU E 

RE  A  C  bEs*0,NCA3€»NDP,GAMA 
!♦  CNL  AAE.Lf  ,0)  GO  TO  -SSJO 
°  R  I  \'  T  b a  J 

0 ;.  ?  <1 ,  N  C  *  S  t  ,  NOP  ,  G  t  »  A 

'  e  >o 

t  c  n « 3 , 

6S J  X ( I ) *0. 

SET4«GAMi/AVG*  ♦  1  ,  I  UMAbCAMA* Aye* 

C  ---  «EA3  Ik  *lPHA,TENO«  BEGINING  AND  END  OF  TI“E  CALCULATION 

OESEECTIVELT 

C  DT«  f !  Hf  INTERVAL  OF  PRINT  OUT 

C  VC"I'<P*CT  VELOCITY 

EMIIJbTJME  SI7ES  OF  INTEGRATION  IN  RUN6E»KUTT A  mFTMQO 


3? 


BEST  AVAILABLE  COPY 


,-r 


.r 


'Y* 


READ  5050  ,  ALPHA, TEND. DT,  V0 
PRINT  5050,  ALPHA, TEND, DT,  V3 
NTTa 1  S  ALPHA0B ALPHA 
READ  5050,  CEM(I),Ial,6) 

PRINT  5050,  (EM(I),I*1,6) 

c  ---  read  in  impulse  shape 

SF0 ( I ) *SF ( I J 

C  K*NO «  OF  TIME  STATIONS  OF  THE  PULSE  SHAPE 

C  TM(J)=  TIME  FROM  IMPACT  BEGINS  OF  STATION  J 

C  CP(J)*RATIO  OF  FORCE  AT  STATION  J  TO  THE 

CC  INPUT  FORCE 

READ  5060,  K , ( TM ( J ) , CP  ( J ) , J*  J  ,  K ) 

PRINT  5060,  rf,(TM(J),CP(J),Jsl»K> 

IF(NDP.LE.0)  GO  TO  625  S  DO  622  1*1, NDP 
C  —  READ  IN  SPRING  CONSTANT  DK  AND  DAMPING  RATIO  DPC 
C  FOR  REPRESENTING  RUBBER  MATERIAL  PROPERTIES 

C  Jl  =  JOINT  NUMBER  OF  THE  ASSOCIATED  MASS 

C  J2*  JOINT  NUMBER  OF  THE  SUPPORT, OR 

C  RIGID  SUPPORT  IF  J2*0 


READ  6200, Jl, J2,0PC,DK 

PRINT  62C0,Ji,J2,DPC,DK 

JlsJi*3-2  $  JP1(I)*MACH(NQ,IR,J1) 

S  J*JP1(I) 

JP2(I)=3 

IF(J2.LE,0  )  GO  TO  621 

J2*J203-2  $  JP2(I)*MACH(N0,IR,J25 

S  DK=5D ( J) 

621 

DP(I)*DPC*2.*SQRT(DK/TR(J)) 

622 

SK  ( I  ) aDK/TR ( J) 

CONTINUE 

625 

CONTINUE 

696 

JNT3*  JNT  +  3  S  JNT 6* JNT *6 

HbEMCNTT) 

612 

DO  612  1*1, IH  S  MO(I)»i 
SF(I)=SF0(I) 

700 

MTF  =  0 

MNFa  1 

00  700  1*1 , JNT6 

VJ£I)*0, 

710 

DO  710  I*2,JNT6,6 

YJ(I)=V0 

KSNQ(I)  $  i4  1 

DO  720  1*1  TNT 3 

IFCI.NE.K)  GO  TC  /20 

12*11*2-1  S  3*12+1  $  K2*K*2-' 

$  K3*K2+1 

YJ(I? *YJ£K2>  S  YJ(I3)=YJ(K3!  S 

11*11+1 

720 

IF(Il.GT.IR)  rJ  TO  725  S  K=NQCI1) 
CONTINUE 

725 

CONTINUE 

730 

LN1*0 

OMEGA* ALPHA^DT 

IF(nMEGAsLE.TEND)  GO  TO  727 
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??e  •••  :•  *■ 

r27  > 

PRINT  fc’04-  Okiil,A 

r'.lL  RG  T  UIPHA, OMEGA, H,YJ, OF) 

:j  7«0  :=2,n.p  s  tisi/2  i  %  nsttodn 

O F G s C r  '  ■  D/J34..-: 
lF£AB5C0*r!;),LF.X(Ii;  '  GO  70  73 1 
*(Il,«DFG  S  V(H)sOMEGA 
736  CONTINUE 

:i-Cl3+?3/3  S  I3»!3-Il*3+3 
14*1/2  i  I5»IAI*1 
IF(I3.E0.1)  GO  TO  731 
GO  TO  (710.733,7355 

731  MNC  =  1  $  IF c Y J Cl  ).LT,T,e  .AND,MG(I4) .tt, 1 )  GO  TO  740 

WTF  =  1 

IF(MQ(IU) ,EQ.l)  732,73a 
73?  Mi2iI4)=2  S  MNTsZ 

733  5F£l4)«a,  S  GO  TO  7 a0 

73a  IF(YJ(I  ),GE,0.)  GO  TO  732 
MG(I4)=3  $  MNF=3 

735  SF(I4)s3F0(Ifl) 

740  PHINT  6100,11,13  , Y J ( 1 2 )  ,  V J C  I) , DF ( 1 , 1)  ,DFG 

IF C AB9(YJ(1) ).GT,100.)  GO  TO  750 
ALPHA=OMEGA  s  GO  TO 
750  IF ( NTT  ,GE,M  GO  TO  600 
ALPHAnALPHAe 
NTT3NTT  +  1  5  GO  TO  696 

5010  FORMAT ( 8F 1 0 , 0  3 
5020  FORMAT( 1 3 , F7 . 2 , 7F  10 . 2 ) 

5030  FORMAT (215,  61 1 , 14, 6F 10,3) 

5040  FORMAT ( 1615) 

5042  FORMAT (4F10.3,/,OF10,3) 

5050  FORMAT (8F10.5) 

5060  FORMAT (1 10, 10X,6r 10.3,/, (8F10a3)) 

5070  FORMAT  (20X,6Fl'?,4) 

5080  FORMAT(I10,4E12.3) 

6001  FORMAT ( 1H1 ) 

6304  FORMAT {  5X,*  T**,F10,5) 

6103  FORMAT (2I5,3E12,3,  2F9,2) 

6200  FORMAT(2I5,3F10,3) 

7000  F0RMAT(12F11,2) 

800  STOP  S  END 


SUBROUTINE  RGKT(A1,A2,A3,YJ,D) 
c  ***  runge-kutta  method  *** 

COMMON  N 

DIMENSION  0(5,80) ,YJ(S0),W( 10), A (10),PHI  (80)  ,8(4,3), 

1  X  ( 10) » Y  (5»8l4> 

C  .  COEFFICIENTS 

W'(l)  =  tN(«)»i./6.  S  W(2)»h(;)=4i/3.  4  kk»4 

A(2)«A(3)»,5  S  A  ( 4 )  *  1  « 

8(2,1) =8 (3, 2)*. 5  S  B(3,l)nB( 4, l)3B{4,2)a0,  S  B(4,3)=l 
IQ  s  4  S  FCT  ■  19./270, 

ALPHA  =  Ai  $  OMEGA  *  A2  $  H  *  A3 
IQbKK 

I QM 1  8  IJ»1  $  IQP1  s  IQ+t 
ISTP  30  $  SIGN  b  l, 

IF (  H.LT.0,  )  SIGN  r  -1, 

X(l)  3  alpha 
DO  3  1=1, N 

3  Y  ( 1 , 1 )  3  Y  J  (  I ) 

4  MM  =  1  $  IFLG  *  0 

h  M  b  MM  %  MM  3  M  ♦  1 

I F (  MM.GT.IQPl  )  MM  a  1 
X (MM )  a  X ( M )  ♦  H 
TEST  a  OMEGA  -  X(MH) 

TEST  1  =  TEST/OMEGA 
IF  QUOTIENT  OVERFLOW  7,8 

7  TEST  1  =  TEST 

8  IF (  A0S(TEST1).LT.1,0E-10  )  GO  TO  12 

I F (  SIGN*TE3T1  )  9,12,13 

9  TEST2  =  OMEGA  -  X(M) 

11  H  a  JEST2  S  IFLG  ■  0 
X (MM)  s  X(M)  ♦  H 

12  ISTP  *  1 

13  XJ  3  X ( M ) 

00  14  I  a  1 , N 

14  YJCI)  a  Y  (M, I ) 

C  .......  RUNGE-KUTTA  PROCEDURE 

DO  25  K*1 , KK 

IF (  K.EQ.l  )  GO  TO  22 

XJ  3  X (M)  ♦  H* A (K ) 

DO  15  1=1 ,N 

15  PHI  Cl )  3  0, 

KM1  s  K  -  1 
DO  19  J* 1 , KM  1 

19  P'il(I)  a  PHI  ( I )  ♦  H*B(K,J)*D(J,I) 

20  Y  J ( I )  k  Y (M, I )  ♦  PHI ( I ) 

22  CALL  OERFCN(XJ,YJ,K,D) 

25  CONTINUE 

C  ........  RKAMSUB  ENTRY  POINT 
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c  - - -  RKAMSU8  EXIT  POINT 

A  2  =  x  ( y  v  3  ?  i  3  =  *- 

00  1 0 5  I s 1 , N 

IMS  Y  J  ( I )  =  YCMM,I) 

501*  RETURN 
E  n  0 

CALL  Dt.RFCN(Al,YJ,l,Oi 
00  27  1  =  1, N 
27  PHICI)  =  0. 

DO  30  I  - 1 ,  N 
DO  29  K*lfKK 

29  PHICI)  =  PHI  ( I )  t  H*WCK)*0(K,I) 

30  Y(MM,I)  =  Y  CM , I  )  ♦  PHICI) 

IF  C  ISTP.EQ.l  )  GO  TO  100 
30  TO  6 

100  CONTINUE 

SUBROUTINE  DERFCN  C  X,? ,  YJ,  M,  OF) 

C  ***  DIFFERENTIAL  EQUATIONS  *** 

COMMON  N,S(120,120),SMKC40,40) 

COMMON  TMC20)  ,CPC20),SFC120),HOC120)fLM,LN2,SLOP 
COMMON  NDP,  JP1C50),  JP2C50)  ,DPC56) 

COMMON  GAMA, SK, BETA  , MTF 

DIMENSION  SQC40),OFC5,80),YJC80),SUC40),3KC120) 

Nl«N/2  $  1FCLN1.EQ.0)  GO  TO  10 
IFCXJ.LE.TMCLNl) .ANO.MTF.EQ, 0)  GO  TO  40 
IFCXJ.LE.TMCLND)  GO  TO  15 
10  lnulnui 

15  CONTINUE 

DO  20  1*1, N1 

20  SQ(I)sHGCI)+CPCLNl)*SFCI) 

DO  30  1*1 , N 1  S  SUCI)»0.  S  DO  30  J*1,N1 
30  SUCI)aSUCI)vAIJCS,I,J,2)*SQCJ) 

40  CONTINUE 

DO  60  1*2, N, 2  %  K*I-1  S  OF C M, K ) *Y J C I )  S  11*1/2 
SUM1*0, 

DO  50  J*  1 , N1  S  J1=J*2-1 

50  SUMlaSUMU8MKCIl,J)*YJ(  Jl)  *BET  A 

sum;=sumugama*yjci) 

DF(M,I)*SUCI1)-SUM1 
60  CONTINUE 

IFCNDP.LE.0)  GO  TO  60  S  DO  70  Iot,N[)P  $  I1*JF1CI)*2 

DVcYjCIl  ) *DP ( I )  $  I2aJP2CI)*2 

IFCI2.LE.0)  GO  TO  65 

OV*(YJCI1)-YJCI2))*OPCI)  S  0FCM,I2)*0FCM,I2)40V 
GO  TO  70 

65  I3«1 1-1  S  0V«DV*SK(1)*YJCI3) 

70  DFCM,I1)=DFCM,I1)-DV 

80  CONTINUE 

MTFag 

RETURN  S  END 
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SUBROUTINE  RGKT(Ai,A2,A3,YJ,D) 

C  ***  RUNGE-KUTTA  METHOD  *** 

COMMON  N 

DIMENSION  D(5,80),  YJ(80),W(10)  M  (10)  ,PHH80),0(4,S} , 

1  X  C 1 0  J  r  Y (5,80) 

C  ------  COEFFICIENTS 

W(l)s*(4)»l,/6.  S  W(2)sW{S)  =  1./3,  S  KK*4 

A(2)sA<3)»,5  $  A  (  4 )  *  1  , 

6(2, l)c@(3, 2)*. 5  S  B (3, 1 ) s0 (4, 1 ) aB ( 4, 2) s0#  $  B(4,3)e! 

IQ  =  4  S  FCT  ■  19./270. 

ALPHA  =  A1  S  OMEGA  *  42  $  H  z  A3 

ICcKK 

IQM1  «  IQ-1  S  IQP1  =  IQ+1 
ISTP  =0  S  SIGN  c  i, 

IF  C  H.LT.0.  )  SIGN  *  -1, 

X(i)  s  ALPHA 
DO  3  1*1, N 

3  Ytl, I)  *  Y J ( I ) 

4  MM  s  1  S  IFLG  *  0 

6  H  MM  %  MM  ■  M  «-  1 

IF C  MM.GT.IQP1  )  MM  =  1 
X  CMH  >  »  X  (M)  ♦  H 

TEST  =  OMECA  -  X(MM) 

TESTi  e  TEST/OMEGA 
IF  QUOTIENT  OVERFLOW  7,8 

7  TESTI  e  TEST 

3  IF (  A0S(TE8Tt),lT,l,0E*10  )  GCl  TO  12 
IF(  oIGN^TESTl  )  9,12,13 

9  TEST2  =  QMFGA  -  V(M) 

It  H  a  TEST2  S  XFLG  a  0 
X (MM)  *  X ( M )  ♦  H 

12  ISTP  ■  1 

13  XJ  s  X ( M ) 

00  14  I  a  1 , N 

14  YJ(I)  e  f  (  M  ,  I  ) 

C  .......  RUNGE-KUTTA  PROCEDURE 

DO  25  K*1 , KK 

IF (  K.EQ.l  )  GO  TO  22 

XJ  a  X(M)  +  H*A(K) 

DO  15  1*1, N 

15  PHIU)  s  0, 

KM1  s  K  -  1 
DO  19  J»1,KM1 

19  PHI(I)  «  PHI  (I)  ♦  H*B(K,J)*D(J,I) 

20  Y J ( I )  x  Y (M, I )  ♦  PH  I ( I ) 

22  CALL  DERFCN(XJ,YJ,K,D) 

25  CONTINUE 

C  .......  RKAMSUB  ENTRY  POINT 
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1 

i  SUBROUTINE  CONF  CA„ “R .MR,ID,LMS,T) 

DIMENSION  Ad?*'. ,  MKd20),T(120) 

DC  10  1*1,  ID  %  1 1 5' M R  (  !  ) 

10  a  ( 1 1 )  *  a  ( 1 1 )  +  a  ( i.  ■  i *  *  t  c ; ) 

AlLMS)s0.0 
RETURN  S  END 

SUBROUTINE  TPAn-5'--  ,  MR  ,  MC  ) 

DIMENSION  AC.j.a  )  -is  i,(0 
DO  20  1*1, NR  a  O'-  .0  J»'1,NC 

20  6(J,l)*A(X,Jj 

RETURN  s  END 

SUBROUTINE  atG>'  iNfl.  Trt,*7#MHI) 

0IMEN3I0N  NOU20)  ,MV(1<‘0),NEW(120) 

Ii*0 

DO  20  1*1, IR  ^  11*11*1 

'  DO  10  K* 1 » MN  l  IF(I*EO.MV(K))  15,10 

15  11*11-1  S  GO  TO  20 

10  CONTINUE 

NEW(I1)*NQ(IV 

r.  20  CONTINUE 

IRbIR-MN  s  do  30  1*1,  IR 
30  NO ( I )  =NEW (15 

RETURN  S  ENO 

SUBROUTINE  SYMINV(A,M) 

DIMENSION  AO 20  - 120)  , T ( 120) 

MisM-i 

DO  12  Irl, Ml  S  TCI )rg.  4  DO  i«  J*i,Hi 
IFCI.GE.J)  2, 4 
2  B* A ( I , J )  S  GO  TO  10 

R  B«A ( J , 1 5 

10  T ( I ) *T  C I ) +B*A (M, J) 

A2*0.  S  DO  20  1*1, Hi 
20  A2«A2+TC15*ACM,I5 

A2*A(M,M)-A2  S  A2*l,/A2  S  A(M,M)*A2 
DO  30  1=1, Ml 
30  A(M,I)*-TCI)*A2 

DO  40  1*1, Ml  S  00  40  J*I , Ml 

40  A(J,I)sA(J,I)-A(M,Z)*T(J) 

RETURN  *  END 
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SUBROUTINE  REMOF(A,IR,MV,MN) 

DIMENSION  A( 122) ,MV( 120) 

11  =  2 

DO  22  1  =  1, IR  $  11=11+1 

DO  12  *=1,MN 

I  *r  f  !  »  E  3  ,  M  V  i  K  ;  i  1 1 

I I  s  1 1  -  t  $  ■"  o  •  C  0  ‘ 

CONTINUE  S  ) 

CONTINUE 

RETURN  S  END 

SUBROUTINE  MMUI.  '*,B,C,NR,N12,NC) 

DIMENSION  a  C  6 ,  e ) ,  3  ( 6,  6 ) ,  C  ( h  ,  t> ) 

DO  2U  1  =  1, NR  S  DO  20  J«1,NC  S  CCI,J)*0. 
DO  20  K=1,N12 

CCI, J)=C(I,J)+A(I,K)*B(K,J) 

RETURN  5  END 


SUBROUTINE  ENDFOC  ( EF , NGO , OF, SL ) 

DIMENSION  EFC6)  ,DF\*to) 

EF ( 2 ) =0 ,  S  EFv5)*0. 

GO  TO  (20,10)  t NGO 
B=SL-DF (2) 

EF(l)  =  CF(n*B  **2*(U  +  2,*DF(2)  )/Sl**3  S  EF ( 4) oDF (P)-EF ( 1 ) 
w?  =  0F(l)*DF(2)*B/Sl**2  %  F.F(3)=-B*wP  S  EF  (b )  *DF  (  2 )  *WP 

GO  TO  30 

IF(OF(3),EQ.0.J  0F(3)*8L 

A  =  DF ( 2 )  $  B  =  SL-DF(3)  S  WP=  DF ( 1 ) / (1 2 ,  *SL*Sl ) 

EF(3)*«P*,,(SL-A)*'*3*(SL  +  3.*A)-B**3*(«.*SL-3.*B)) 

EF(6)s-WP*( (SL-6)**i*(9L  +  3.*B)-A**3*(4.*SL-3. *A) ) 

SHbDF (31-DF (2)  S  +  SH '2 . ) +DF ( 1 ) 

EF(1)s(WM+CF(3)+EF(6))/SL  8  FF(<|)=SH*DF  ( 1  )-EF(l) 

FF(3;=-EF(3)  J  EF (&)*»EF (6) 

RETURN  S  END 


SUBROUTINE  ADDMS ( W , M , S , SD , IR , I J, IK ) 
DIMENSION  W(6,6) ,S(120, 120) ,SOC120) , I JC120) 
DO  5  1=1 , IR 
S ( I » I ) =SD  (I ) 

DO  30  I = 1 , M  S  I1«IJ(I) 

DO  20  J « I , M  $  J1  =  IJ(J? 


IF(Jl.GE,IltANO,IKtEO.l) 
IF ( Jl.LE,IltAND.IK,FQ,2) 
GO  TO  15 

S(Il,Jl)aS(Ii,Jl)+W(I,J) 

s(ji,ii)*s(ji,n)+wa,j) 

CONTINUE 

SDt:i)*s(ii,in 

CONTINUE 

RETURN  $  END 


GO  TO  10 
GO  TO  10 

S  GO  TO  20 
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SUBROUTINE  MUL1 (A,B,C, IR) 

0 1  MENS  ION  A(6,6),B(6),C(6) 

DO  10  1=1, IR  $  C(I)*0, 

DO  10  J*1,IR 

10  C(I)*C(I)+A(I, J)*8(J) 

RETURN  $  END 

FUNCTION  MAC.H(NQ,IR,KI) 

DIMENSION  NQ ( 1 20  D 
DO  iO  1  =  1,  IR 

!F(NQ(I).Nfc,Kl)  GO  TO  10 
MACH=I  S  GO  TO  20 
10  CONTINUE 

20  CONTINUE 

RETURN  $  END 

SUBROUTINE  C0NMS6(S,M,MR,ID,LM5,TR) 

DIMENSION  8(6,6) ,MRC 120) ,TR( 120) ,T1 (120) »T2( 120) 
IK*1 

20  00  30  1=1, M 

IFCI. GE.LMS, AND, IK. EQ.l)  21,22 

21  T 1 ( I ' =S (LMS, I )  S  GO  TO  30 

22  T1(Ij=S(I,LM8) 

30  CONTINUE 

DO  110  1  =  1, ID  S  1 1 «MR  Cl) 

DO  40  J=1,M 

IFCJ. GE.il, AND, IK, EQ.l)  31,32 

31  A=S(I1,J)  S  GO  TO  <10 

32  A=S(J,I1) 

40  T2(J)rA>TR(I)*T1(J) 

DO  50  J  =  1 , ID  S  J 1 =MR ( J) 

50  T2CJ1)=T2(J1)+T2(LMS)*TR(J) 

L 1  *  I  * 

DO  100  K  =  1,M  S  IF (LI  ,LE,0)  GO  TO  70 
DO  60  L»1»L1 

IF(K.EO,MR(L)J  GO  TO  100 
60  CONTINUE 

70  IF(K,GE.Il.ANO,IK.EQ.n  80,90 

80  S ( 1 1 , K ) =T2 (K)  $  GO  TO  100 

90  S(K,I1)=T2(K) 

100  CONTINUE 

J.10  continue 

DO  150  1*1, M 

IF(I,GE.LM8,AND,IK,EQ,l)  120,140 
120  SCLMS,n=0,  I  GO  TO  150 

140  5(1  , LMSi «0 , 0 

150  CONTINUE 

DO  160  1*1, M  S  DO  160  J»I,M 
163  S(J,I)=S(I,J) 

RETURN  S  END 


SUBROUTINE  C0NM3(S,SD,M,mR, ID,LMS,TR,  IK) 

DIMENSION  S(120,120),SD(120),MR(120),  TR ( 1 20) , T1 ( 1 20) , 
1  T2(120) 

6  DO  7  1=1, M 

7  Sd.I)aSD(l) 

20  DO  30  1*1, M 

T1 C • )=ATJ(S 
30  CONTINUE 

DO  110  1*1, ID  S  IlsMRU) 

DO  40  J*1,M 
A*AIJ(S,I1,J,IK) 

T2( J)*a*7R(I)*T1 (J) 

40  CONTINUE 

DO  50  J= 1 , I D  S  Ji*MR(J) 

50  T2CJ1 )*T2CJ1)TT2(LMS)*TR(J) 

L1«I-1 

DO  100  K* 1 , M  S  IF (L  1 , LE . 0 )  GO  TO  70 
DO  60  1*1, LI 

IF(K.EO.MRCL)  )  GO  TO  100 
60  CONTINUE 

70  IF(K.GF.I1.AND,IK#EQ,1)  GO  TO  80 

IF ( K . LE  « I  1  #  AND , IKfEO,2)  GO  TO  80 

GO  TO  90 

80  S(I1»K)»T2(K)  *  GO  TO  100 

90  S(K,I1)*T2(K) 

100  CONTINUE 

110  CONTINUE 

SaMS,LMS)*C. 

DO  150  1*1, M  S  SD(I)«S(I,I) 

150  CONTINUE 

return  *  END 

SUBROUTINE  REMOMS(S,SD,  IR, MV , MN , IK ) 

DIMENSION  5(120, 120),SO( 120) ,MV( 120) ,T( 120) 

DO  2  1*1, IR 
2  S(I,I)=SO(I) 

11*0 

DO  50  1*1, IR  S  I1*I1M 

DO  10  K*1,HN  $  IF(I.NE.MV(K))  GO  TO  10 
1 1  *  1 1  ••  1  $  GO  TO  50 

10  CONTINUE 

K*0 

DO  30  J*I,IR  J  K*K+1 

DO  20  L*1,MN  $  IF (J.NE.NV(L) )  GO  TO  20 
K*K-1  $  GO  TO  30 

20  CONTINUE 

T(K)cA,lJ(SrX,Jf  IK) 

30  CONTINUE 

DO  40  J* 1 , K  5  L*I1+J-1 

IFU..GE,  Il.AND.lK.EQ.l)  GO  TO  32 
IF(L»LEtIl«AND.IK,EGI»2)  GO  TO  32 
GO  TO  23 
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GO  TO  40 


32  S(I1,L)«T(J)  S 

33  S(L,I1)«TCJ) 

40  CONTINUE 

8ocii)«s(ii#ir.) 

50  CONTINUE 

RETURN!  S  END 

FUNCTION  AIJ(S,I1,I2,I«) 

DIMENSION  S(120,i2e) 
IF(I2.GE.Il.AND.IK.EG.n  GO  TO  1 
IF(I2.LE,IltAND.IK.EQ.2>  GO  TO  1 
GO  TO  2 

1  AIJ»S(I1,I2)  i  GO  TO  3 

?  UJ3SCI2,I!) 

3  RETURN  S  END 


SUBROUTINE  MIV(A,U,NM) 

DIMENSION  A(6,6),U (6,6) 

00  9001  I»1,NM 
00  9001  J=1,NM 
U(I, J)*0, 

IF  (i.tQ.J)  U(I,J)=1.0 
9001  CONTINUE 

EPS30. 0300001 
DO  9015  1*1, NM 
K3I 

IF  C 1 -NM  J  9021  r  9007  »  90  c 1 
9021  IF  C A  €  I » I J-FPSJ  9005,9006,9007 

9005  IF(-A(I,I)-EPS)  9006,9006,9007 

9006  ‘<ak'+l 

DO  9023  ^ 3 1 , NM 
U(I,J)sUCI,J)+U(K,J) 

9023  ACI,J)  =  A(I,J)4.A(K,J) 

GO  TO  90tl 

9007  DIV3A(I,I) 

DO  9009  J*1,NM 
UCI, J)«U(I, JJ/OIV 

9009  4(I,JisA(I,J)/DIV 
DO  9015  MM* 1 , NM 
DElT<*A(MM,I) 

IF  (ABS(DELT)-EPS)  9015,9015,9016 
9016  IF  C MM- 1 5  9010,9015,9010 

9010  DO  9011  J*1,NM 

U(MM, J)«U(MM, J)-U(I# J)*DELT 

9011  A (MM, J)SA(MM, JD-ACI, J)*D£LT 
9015  CONTINUE 

RETURN 

END 
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APPENDIX  B 


Cushioning  System  Design 
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Cus hjoning  System  Desig i: i 
The  weight  distribution  of  an  M~  :7  3/4  ton  truck 
with  1500  lb.  simulated  load  of  sandbags  and  the  arrange¬ 
ment  of  cushioning  forces  are  shown  in  F:.g.  B.l.  In  order 
to  have  uniformity  in  the  crushing  of  the;  honeycomb  cushion¬ 
ing,  the  centroid  of  the  honeycomb  area  or  the  center  of 
cushioning  forces  should  be  located  at  the  center  of  gravity 
of  the  loaded  vehicle.  Wheels,  differencials,  gear  reducer 
and  transmission  can  be  cushioned  independently  and  the 
cushioning  forces  for  those  masses  are  calculated  as: 

(Design  acceleration  level  17.5  g) 

Wheels  and  tires  (each 

Differentials  (each) 


Gear  Reducer 
Transmission 

The  remainder  of  the  structure  is 
throucf’a  the  truck  frame. 


F  =  18.  F 
w 

•  350  = 

6480  lb. 

F  =  F 

fd  rd 

=  18.5 

x  480 

=  8890  lb 

F  =18.5 

gr 

FTr  =  18.5 

x  300  = 

5550  lb 

x  200  = 

3700  lb 

cushioned  by 

Pl'  2F2 

and  2F4 

or 


EF  =  0  (Balance  of  cushioned  forces  with  mertia.1 
forces) 

Fx  +  2F2  +  2F4  =  (1500  +  1060  +  570  +  600  +  300)  x  18 

F1  +  2F  2  +  2F4  =  74'500  (B.l) 

EM  =  0  (Sum  of  moments  of  cushion  forces  about  CG) 

2F .  x  80.5  +  (2F  f  F  ,)  X  50.5 
4  w  rd' 


=  2F2  +  FTr  x  31  +  (2Fv;  +  Ffd)  x  63 


or 


+  81  x  F,  +  5  x  F 

1  gr 


161F4  -  2F9  -  81F1  415,325  (B.2) 

Since  we  have  three  unknowns,  F^,  F2,  and  F^  we  need 
one  more  constraint  equation.  This  can  oe  achieved4by  intro¬ 
ducing  an  artificial  hinge  at  the  internal  cushioned  point  F2< 
If  we  take  the  moment  about  this  point,  we  get: 


2F4  x  81.5  =  (1500  +  1060)  x  45  x  18.5  (B.3) 


Thus  F]_  =  20,600  lb. 

F2  =  13,900  lb. 

F4  =  13,100  lb. 

The  sizes  of  the  pads  required  for  the  cushioning  forces 
calculated  by: 


are 


where 


A  =  required  area  of  honeycomb  under  'cushion  force  F. 
S  =  average  stress  of  honeycomb  (6400  psf) . 

From  a  work  energy  balance,  the  stack  height  z  required  to  . 
provide  the  volume  of  honeycomb  necessary  to  cushion  the  vehi¬ 
cle  is  determined  from: 


where 


H  _  10  x  12 

Z  Ge  17.5  x  0.7 


9.7" 


H  =  equivalent  free  drop  height  (10  ft.) 
G  =  design  acceleration  level  (17.5) 
e  =  design  strain 


) ■» 
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